Self-consistent symmetries in the proton-neutron Hartree-Fock-Bogoliubov approach 



o 
o 

(N 

o 
O 

00 
(N 



> 
o 

in 
d 

O 



S.G. Rohozinski/ J. Dobaczewski,^' ^ and W. Nazarewicz^^^'"* 

'institute of Theoretical Physics, University of Warsaw, ul. Hoza 69, PL-00681, Warsaw, Poland 
'^Department of Physics, University of Jyvaskyla, 
P.O. Box 35 (YFL), FI-4OOI4 University of Jyvaskyla, Finland 
'^Department of Physics and Astronomy, The University of Tennessee, Knoxville, Tennessee 37996 
^Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, Tennessee 37831 

(Dated: October 28, 2009) 

Symmetry properties of densities and mean fields appearing in the nuclear Density Functional 
Theory with pairing are studied. We consider energy functionals that depend only on local densities 
and their derivatives. The most important self-consistent symmetries are discussed: spherical, axial, 
space-inversion, and mirror symmetries. In each case, the consequences of breaking or conserving 
the time-reversal and/or proton-neutron symmetries are discussed and summarized in a tabulated 
form, useful in practical applications. 

PACS numbers: 21.60.Jz, 21.30.Fe, 21.10.Hw, ll.30.Er, 71. 15. Mb 



I. INTRODUCTION 

The nuclear Density Functional Theory (DFT) [E[l,[| 
is a theoretical tool of choice for describing complex, 
open-shell nuclei, for which the dimension of the config- 
uration space becomes intractable for other methods of 
theoretical nuclear structure, such as ab initio and con- 
figuration interaction (shell model) techniques [1]. The 
main building blocks of nuclear DFT are the effective 
fields, often represented by local proton and neutron den- 
sities and their derivatives. When compared to the elec- 
tronic DFT, the unique features of the nuclear variant 
are: (i) the presence of two kinds of fermions, protons 
and neutrons; (ii) the essential role of nucleonic pairing; 
and (iii) the absence of external potential and the need 
for symmetry restoration in a self-bound system. 

At the heart of the DFT lies the energy density func- 
tional H that is built from the nucleonic intrinsic den- 
sity matrices. The requirement that the total energy 
be minimal under a variation of the densities leads to 
the Hartree-Fock-Bogoliubov (HFB; or Bogoliubov-de 
Gennes) equations. They form a set of non-linear integro- 
differential equations which has to be solved iteratively 
for the self-consistent densities. 

The quasi-particle vacuum associated with the DFT 
solution is a highly correlated state. This is partly be- 
cause the DFT description is performed in a frame of 
reference attached to the nucleus, the intrinsic frame, in 
which the nuclear mean field may spontaneously break 
the symmetry of the original Hamiltonian, or energy den- 
sity. While the resulting deformed solutions do not obey 
symmetries present in the laboratory system, they ac- 
quire lower binding through long-range polarization ef- 
fects. Additional correlations may be gained by means 
of symmetry restoration. Such a strategy, rooted in the 
nuclear Jahn- Teller effect Q, has proven to be very 
effective in nuclear mean- field calculations 0, Q ■ 

Since the symmetry breaking is essential for deter- 
mining the optimal mean field of the nucleus, the self- 
consistent symmetries present in the model [l^, [HI [H, 



Il3|,ll4|,ll5| may often determine physics. A self-consistent 
symmetry (SCS) is a unitary or antiunitary transforma- 
tion iS, which commutes with the HFB Hamiltonian. Due 
to self-consistency, S also commutes with DFT densities. 
The associated DFT energy density H is referred to as 
symmetry invariant [l^ : 



(1) 



The above definition of SCS that can often be found in 
the literature Q is too limiting when it comes to the 
DFT. Indeed, invariance of the energy density itself is 
not a prerequisite for the invariance of the energy density 
functional (EDF). Actually, the EDF is invariant with 
respect to S also when the energy density is covariant 
with S, i.e.. 



H^ir) =n{S+rS). 



(2) 



The energy density that meets ([2]) is referred to as 
symmetry covariant, see discussion in Appendix A of 
Ref. [I61]. The existence of SCS has a profound impact 
on self-consistent solutions. If the initial density matrix 
employed at the first iteration of HFB equations contains 
a SCS, then that symmetry will propagate through to the 
final DFT solution. Therefore, the introduction of SCS 
restricts the generality of the self-consistent density ma- 
trix and may lead to an erroneous estimate of the DFT 
energy and deformation of the system. 

A considerable literature exists on nuclear collective 
modes associated with spontaneously broken symmetries. 
Ground-state deformations of nuclei, including those in 
the pairing channel, have been reviewed in Refs. [l7l.lT8l|. 
In the presence of angular momentum, new deformations 
of magnetic character may appear in the rotating nuclear 
mean field. High-spin particle-hole (p-h) and particle- 
particle (p-p or pairing) deformations, both isoscalar and 
isovector, have been discussed in detail in Refs. (Tol. |20|. 
which also contain a general discussion of spontaneous 
symmetry-breaking phenomena in rotating nuclei (see 
also for a recent update). 
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Symmetry properties of Hartree-Fock (HF) densities 
were studied in Ref. [22] in the context of a double 
point group D™ that contains three mutually perpendic- 
ular symmetry axes of the second order, space inversion, 
and time reversal. The associated symmetry-breaking 
schemes have been outlined in the following Ref. [23| . 

In this paper, we extend the discussion of Refs. [2^[23j 
to the pairing channel and transition densities using the 
coordinate-space HFB theory which incorporates an ar- 
bitrary mixing between protons and neutrons in the p- 
h and p-p channels The constraints imposed on 

DFT densities by the time-reversal and proton-neutron 
(p-n) symmetries are studied for various spatial geome- 
tries (spherical, axial, space-inversion, mirror, and D2h)- 
The details pertaining to the p-n HFB theory in the Local 
Density Approximation (LDA) can be found in Ref. [13] . 
Throughout our work, we refer to this previous paper as 
I and to, e.g., Eq. (1) therein as (I-l). However, to make 
the present paper self-contained, some of the definitions 
given in I are also repeated here. 

Our interest in a general HFB formalism that incor- 
porates an arbitrary mixing between proton and neu- 
tron quasiparticles is motivated by numerous phenom- 
ena that are present in medium-mass and heavy nuclei 
with N^iZ . These include: p-n pairing correlations; al- 
pha decay and alpha clustering; local increase in binding 
(Wigner energy); interplay between isospin T=0 and 1 
states in N=Z nuclei at low and high angular momenta; 
isospin mixing and mirror symmetry breaking; and beta 
decays and superallowed beta decays in particular, just 
to mention a few. All these cases involve isospin in one 
way or another, and since we are often interested in nu- 
clear states with nonzero spin, proper treatment of time- 
reversal symmetry is crucial. A considerable amount of 
literature exists on a mean-field treatment of Nk,Z nu- 
clei. For instance, for a concise overview of p-n pairing, 
we refer the reader to Ref. ^] . 

The paper is organized as follows. SectionlTTlintroduces 
the local DFT densities and defines transformation rules 
for density matrices. The symmetries of interest, both 
in the position-spin space and in the isotopic space are 
discussed in Sec IHII Symmetry properties of densities 
are studied in Sec. IIVI for (i) spherical symmetry (with 
and without space inversion); (ii) axial symmetry (with 
and without space inversion); and (iii) D2h symmetry. 
When going beyond the mean- field approximation, e.g., 
by using the Generator Coordinate Method or projection 
techniques, multi-reference transition densities appear. 
The associated symmetry properties are summarized in 
Sec. |Vl Previous self-consistent calculations of p-n pair- 
ing are commented upon in Sec. IVII in the context of our 
findings. Section IVIII summarizes the main results of 
our study. Finally, Appendix A is devoted to the gener- 
alized Cayley-Hamilton theorem for irreducible spherical 
tensors and tensor fields. 



II. DENSITIES AND FIELDS 

A. Density matrices and mean fields in tlie p-h and 
p-p channels 

To fix the notation, we begin with a brief recapitula- 
tion of definitions and basic properties of the one-body 
HFB density matrices in p-h and p-p channels, see I for 
details. The p-h and p-p density matrices are defined, 
respectively, as 

p{rst,r's't') = {■^\a+^,^,arst\^), (3) 

P{rst,r's't') = As't'{^\ar'-s'-t'arst\'^), (4) 

where a^^j and a^st create and annihilate, respectively, 
nucleons at point r, spin s=±^, and isospin t—±^, while 
is the HFB independent-quasiparticle state. 
The p-h and p-p density matrices form together the 
projective generalized "breve" density matrix. 



'R.{x,x') 



_ ( p{x,x') 



p{x,x') 



p+ix,x') i-/5^^(x,x') 



(5) 



where we abbreviated the position-spin-isospin variables 
by x={rst} and 1 :— S{x — x') :— S{r — r')Sss'Stt'- 

In I, we found that instead of using the usual anti- 
symmetric pairing tensor it is more convenient to 
introduce the above p-p (or anomalous) density matrix, 
p ([4]). The relation between the standard density matrix 

7?.(x,a:') [3 and ']Z{x,x') is given by a unitary transfor- 
mation I24II 



W = 



1 







(72'f(2) 



(6) 



The quantities expressed in the representation ([6]) are 
indicated by a "breve" symbol in the following. 

Throughout our paper we apply the following nam- 
ing convention. The matrices are denoted by "hat" ; the 
quantities expressed in the representation Q are marked 
by "breve" ; the matrices in a double HFB space are de- 
noted by the calligraphic capital letters. As in I, we label 
space vectors by boldface symbols and their scalar prod- 
ucts by the central dot, e.g., r-V; the components of 
vectors and tensors are labeled with indices a, 6, c, and 
the names of axes are x,y, and z, e.g., r = (rx,ry,rz)- 
Here we note that similarly as in I, the individual vector 
components are also shown in boldface. 

Vectors in isospace (isovectors) are labeled by arrows 
with their scalar products in the isospace denoted by the 
circle, e.g., vow. The components of isovectors are la- 
beled with indices i. A:, and the names of iso-axes are 1,2, 
and 3. Isoscalars are marked with subscript "0", and we 
often combine formulae for isoscalars and isovectors by 
letting the indices run through all the four values, e.g., 
/c=0,l,2,3. 

The symbol ^da:; represents integration over spatial 
coordinates and summation over spin and isospin in- 
dices, and • denotes the matrix multiplication and in- 
tegration/summation "^Ax. The asterisk stands for the 
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complex conjugation. The spin and isospin Pauli ma- 
trices are, respectively, ass' = t<^Is' ^'^la') ^"^^ ^tt' — 
(fjjV , f^f) , fj:f} ) , and the corresponding unity matrices are 

T^°} = 5tt' and ct^^, = 5ss' ■ 

Under TC, the product of time reversal T and charge 
reversal C (a rotation by tt in isospace around the second 
axis), the density matrices ^ and (U become: 



16ss'tt'p*{r 



-U r' 
-U r' 



-t')(7a) 
-t')(7b) 



The symmetries of p and p can be conveniently expressed 
in terms of just the Hermitian conjugation and TC: 



(8aj 
(8b) 



Expressed in terms of spin-isospin components, the 
density matrices can be written as: 



p{rst,r's't') = \pa{r-,r')5ss'5tt' + jSss' pir,r') o Tt^, + \so{r,r') ■ &ss'Stt' + js{r,r') ■ &ss' ° fw , (9a) 
p{rst,r's't') = jpo{r,r')Sss'Stt' + jSss'p{r,r') o t^, + iso(r,r') • &ss'Stt' + js{r,r') ■ &ss' o fw ■ (9b) 



To avoid confusion, the functions of two position vectors 
r and r' appearing on the right-hand sides of Eqs. (O will 
be called the nonlocal density functions or, simply, den- 
sities, unlike the density matrices ^ and ^ appearing 
on the left-hand sides. 

Since the p-h density matrix and the Pauli matrices are 
both Hermitian, according to (1-16) all p-h densities are 
Hermitian as well; hence, their real parts are symmet- 
ric, while the imaginary parts are antisymmetric with 
respect to exchanging r and r' . Similarly, transforma- 
tion properties of Pauli matrices under time-reversal and 
charge conjugation (1-17) make the p-p densities either 
symmetric (scalar-isovector and vector-isoscalar) or an- 
tisymmetric (scalar-isoscalar and vector-isovector) with 
respect to exchanging r and r', see Eq. (1-18). These 
properties are fulfilled independently of any other sym- 
metries conserved by the system; they are consequences 
of definitions of density matrices p and p. 

In the HFB theory with the zero-range Skyrme interac- 
tion [25, 26]j_or in the local density approximation (LDA; 
cf. Refs. [3 [13) J the energy functional depends only on 
local densities, and on local densities built from deriva- 
tives up to a given order, see Refs. [lEIlll for systematic 
constructions. The local densities are denoted by hav- 
ing only one spatial argument to distinguish them from 
the nonlocal densities. Following the standard definitions 
[2^ . [3^ , in the present study we employ definitions of lo- 
cal p-h and p-p densities according to Ref. I. For the sake 
of completeness, we repeat them here: 

• scalar densities: 

- particle and pairing densities: 

Pk{r) = Pk{r,r')r=r', (10a) 
p*(r) = p*(r,r').=w, (10b) 



- p-h and p-p kinetic densities: 

n{r) = [{V-V')pkir,r')]^^^„ (11a) 
f(r) = [(V- V')p*(r,r')],^,„ (Ub) 

• vector densities: 

- p-h and p-p spin (pseudovector) densities: 

Sk{r) = Sk{r,r')r=r', (12a) 
So{r) = SQ{r,r')r=r', (12b) 

- p-h and p-p spin-kinetic (pseudovector) densities: 

T,(r) = [(V- V>fe(r,r')],^,,, (13a) 
fo(r) = [(V- V')so(r,r')],^,„ (13b) 

- p-h and p-p current (vector) densities: 

h{r) = ^^[(V- V')pfe(r,r')],^,„ (14a) 
Ur) = 5l[C^-'^')Po(r-,r')],^,„ (14b) 

- p-h and p-p tensor-kinetic (pseudovector) densities: 

Ffc(r) = i[(V®V' + V'®V)-Sfc(r,r')]^^^„(15a) 
Fo(r) = i[(V®V'-KV'®V)-So(r-,r')],^,,,(15b) 

• tensor densities: 

- p-h and p-p spin-current (pseudotensor) densities: 

Hr) = ^[(V- V')®s4r-,r')],^,„ (16a) 
J>) = li[(V- V')®i(r,0],^,,, (16b) 
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where fc=0,l,2,3, and ® stands for the tensor product of 
vectors in the physical space, e.g., {v(>^ w)ab = VaWi, and 

[{v®w)-z]a=Va{w-z). 

The kinetic, spin-kinetic, and tcnsor-kinctic densities 
are, in fact, equal to contractions of the following 2nd- 
order and 3rd-order tensor densities: 



heir) = [VbVj{r,r')\^^^,, 
fobcdir) = [VbV'Jod{r,r')]^^^,, 



namely. 



f (r 
Tkd{r 
fod{r 

Fkb{r 
Fob{r 



Tkbbir), 

b—x,y,z 

^bbir), 

b—x,y,z 

Y Tkbbdir), 

b—x,y,z 

Y Tobbdir), 



b=x,y^z 

1 

2 



i Y iTkbcc{r)+Tkcbcir)), 

c—x.y,z 

I Y {Tobcc{r)+focbcir)). 



(17a) 
(17b) 
(17c) 
(17d) 

(18a) 
(18b) 
(18c) 
(18d) 

(18e) 
(18f) 



c—x,y,z 



All pseudotensor densities can be decomposed into 
trace, antisymmetric, and symmetric components (1-26)- 

(1-28), that is, into pseudoscalar Jfc(r) and J(r); vector 

Jfe(r) and J(r); and pseudotensor Jfc(r) and J(r) densi- 
ties. 

In the case of the Skyrme effective interaction, as well 
as in the framework of the energy density functional ap- 
proach, the energy functional is a three-dimensional spa- 
tial integral. 



H = / d^rn{r), 



(19) 



of the local energy density H{r) that is supposed to be 
a real, scalar, time-even, and isoscalar function of local 
densities and their derivatives. 

Minimization of the energy functional with respect to 
the p-h and p-p density matrices under auxiliary condi- 
tions 



and 



J d^rpoir) = A 



(20) 



(21) 



leads to the common eigenvalue problem for the gener- 
alized density matrix TZ and the generalized mean-field 
Hamiltonian matrix defined as 



h~X h 



(22) 



with the Lagrange multiplier matrix given by 

A = i(A„+Ap)/+i(A„-Ap)f(3) = Ao/ + A3f(3), (23) 

where A„ and \p are the neutron and proton Fermi levels, 
respectively. 

B. Transformation rules for the density matrices 

A general Hermitian one-body operator in the Fock 
space can be written as 

G = j dV^ / d''rYg{r's't',rst)a+^,,,arsu (24) 



where 



g{r'st',rst) = {r' s't'\g\rst) 



(25) 



is the matrix element of the single-particle operator g act- 
ing in the space of one-body wave functions. Let us now 
consider a unitary transformation U in the Fock space 
generated by G: 



(26) 



where a is a real parameter. By making use of the Baker- 
Campbell-Hausdorff relations, the annihilation and cre- 
ation operators transform under U as: 

U+arstU = J d^r'YHrst,r's't')ar's't', (27a) 

s't' 

U+atstU = /'dV5]{t+(rVi',rsi)a+,,,„(27b) 



where 



u{r's't',rst) = {r's't'\e'"3\rst). 



(28) 



From Eqs. pTj) it follows that the density matrices cal- 
culated for the transformed state U\'^), i.e., the trans- 
formed density matrices, are: 
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fix. 



p^{xi,x\) = As[t'^^dx2dx'2'is'2t'2u{xi,X2)p{x2Tx'2)u{x'-^,x'2)., 

where x = {r, — s, —i}. Using a shorthand notation, Eqs. ([29| can be written as 



+ 



(29a) 
(29b) 

(30) 



where u is defined as 



iL{rst,r' s't') ^ lQss'tt'u*{r - s - t,r' - s' - t') = ^ ^ (-^f,(fj)(~^aJ,,Ou*(rs'Y^ r's'''t''0(zCT^,,,,0(«^^^^^^^^ (31) 



V 



with 



(-if(2))(-za^)it*(lCT!')(if^^J), 



(2)^ 



(32) 



(33) 



cf. Eqs. ([7]). It immediately follows from Eqs. ([30ll that 
the generalized density matrix ([5| transforms under U 



(34) 



where the transformation matrix in the doubled- 
dimension space is defined as 



U = 



ii 
iL 



(35) 



Similar definitions can be introduced for any unitary 
antilinear transformation operator, Uk, which can al- 
ways be presented in the form: 



Uk = UK, 



(36) 



where ?7 is a linear unitary operator and K is the opera- 
tor of complex conjugation in the position-spin-isospin 
representation, in which the basis states \rst) are as- 
sumed to be real: K\rst) = \rst). Let us recall that 
the action of if on a single-particle wave function <i>(rst) 
expressed in a basis (l)i{rst), ^{rst) = ^■Ci0i(rsi), is 
defined as 



K^irst) = $*(rsi) = ^ c*(/)*(rsi). 



(37) 



The generalized density matrix ^ transformed under 
Uk reads: 



where the transformation matrix is 



Uk 



UK 
iiK 



(38) 



(39) 



Uk — uK, 
Uk ~ iiK, 



(40a) 
(40b) 



and K is the matrix complex conjugation operator asso- 
ciated with the position-spin-isospin representation. One 
has to remember that the decomposition (|36[) and defi- 
nition of K (f37|l does not depend on any specific choice 
of the single-particle basis 4)i{rst). The advantage of 
such a choice is that properties of antilinear symme- 
tries (like the time reversal) directly translate into the 
complex-conjugation properties of densities. However, 
other choices of K can be useful when the complex- 
conjugation properties of matrix elements of operators 
in the given basis (f>i{rst) are considered. 



III. SYMMETRIES 

Let us suppose that U (or Uk) is a symmetry transfor- 
mation of the nuclear many-body Hamiltonian, H, i.e., 



UHU^ 



H. 



(41) 



The generalized density matrix (O and mean-field Hamil- 
tonian (|22p . obtained through the minimization proce- 
dure, may, but need not, obey the symmetry U. It can 
be only proved [3] that if (7 is a symmetry of H then the 
transformed generalized mean-field Hamiltonian depends 
functionally on the transformed generalized density ma- 
trix in the same way as the original Hamiltonian on the 
original density: 



(42) 



This means that in order to understand the symme- 
tries of the mean field, it suffices to analyze the sym- 
metries of the underlying density matrix. The nuclear 
Hamiltonian is supposed to conserve numbers of protons 
and neutrons, and to be invariant under space rotations 
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D{a,f3,j), space inversion P, time reversal T, and rota- 
tions in the isotopic space (isorotations) DT{ar,l3r,jT)- 
(Throughout the present paper, to denote rotations we 
use symbols D, d, and V instead of the usual letter R. 
The symbol TZ is reserved for the generalized density ma- 
trix.) In the space-spin-isospin basis, the single-particle 
matrix elements of corresponding operators d, p, t, and 
dr are: 

d'^'^^ir's't^rst) = S{r' - r)5vte'^^'^''^ e'l^^y^"-^ e'"^^^''^ 
^ 6{r' - r)(5t'te^'^'^('')e'^'"('')e'"'^('') 
xas's(<T,Q;/37), (43a) 
p[r's't\ rst) - 6{r' + r)5s's5t'u (43b) 
i{r's't' ,rst) = 6{r' - r){-ialJSt'tk, (43c) 
rfa./3.7. ^r's't', rst) - Sir' - r)Ss'sat'ti^, a^Mr), (43d) 

respectively, where l{r) is the single-particle orbital- 
angular-momentum operator, j{r) = l{r) + ^cr is the 
total single-particle angular-momentum operator, and a, 
/3, and 7 (ofr, /3r, and 7t) are the Euler angles of rota- 
tions in space (isospace). The spin rotation matrix ds's^ 
being the function of the Pauli matrices and Euler angles, 
reads 

a,,,(a,a/37) = (e^^-^e^/^-^^"-^),,, 

(if 7 + ar •• 1 + - z 
= cos - I cos —;^5s's + I sm cr,,. 



Improper rotations supplemented by time reversal 
5c]): they form the group called 0^(3). 



.. P f 1 — a . 7 — 

+ism- ( cos^— +sm^— (T^,, 



(44) 



and the isospin rotation matrix aft (r, Q!r/3r7r) is defined 
analogously. 

Below, rotations by angle tt about the three axes x, 
y, and z, which are called signature operators, will be of 
particular interest: 

f^ir's't', rst) = ^^'"'(r's't' , rst), (45a) 
fy{r's't', rst) = d°''°{r's't',rst), (45b) 
f,{r's't', rst) = d'^°"(r's'i', rst). (45c) 

Products of signature operators and the space-inversion 
operator, which are called simplex operators, correspond 
to reflections with respect to the y-z, z-x, and x-y plains, 
respectively: 

Sx=fx» P, (46a) 
Sy ^ fy • p, (46b) 
Sz=f^»p. (46c) 

Symmetry operations (|43p form well-known group 
structures: 



• Proper rotations (|43ap : they belong to the orthog- 
onal unimodular group in three dimensions SO (3). 



Improper (or mirror) rotations, i.e., rotations (|43ap 
combined with space inversion (|43bp : they belong 
to the full orthogonal group 0(3). 



Together with the group of isorotations (|43dp , the 
symmetry operations (j43p constitute the group 
0^(3)xSO(3). 



• Space-inversion (|43bp together with three signa- 
tures (|45l) and three simplexes (|46)) constitute the 
point group of symmetries of a parallelepiped, 
called D2h, which is of interest for triaxial nuclei. 

Due to the fact that the rotation of a spin-i system by 
27r changes sign of the wave function, one has to, in fact, 
double these groups (see Refs. [13, HI] for details). The 
doubling of groups has no bearing on the results of the 
present study; hence, in what follows we do not refer to 
it. 

Although strong nuclear forces are charge- 
independent, i.e., invariant under rotations in isospace, 
the nuclear Hamiltonian is not, at least because of 
electromagnetic forces. Within the HFB theory, the 
charge independence is additionally broken by the aux- 
iliary condition (|2ip which is manifestly isovector. This 
constraint gives rise to HFB product states which violate 
isospin even for charge-independent Hamiltonians (see 
Ref. 31] for a recent discussion). If the density matrix in 
the p-p channel does not vanish, particle number is also 
violated in the HFB theory. The local energy density 
7i(r) is usually constructed under assumption that it 
should be invariant ([1]) with respect to the time reversal 
T and isorotations Dr, and covariant Q with respect 
to the space symmetries D and P, see Appendix A in 
Ref. [1^. All these symmetries are often spontaneously 
broken in mean-field theories. As discussed in the Intro- 
duction, the problem of symmetries that are conserved 
by H and internally broken by Ti. is, in fact, one of the 
most important elements of a mean-field description of 
many-body systems. 



A. Symmetries in the isotopic space 

As has been discussed in I, the standard case of no 
explicit p-n mixing can be described by the conserved 
p-n symmetry given by 



C3 



-m(f',7r00) = f(3^ 



(47) 



That is, C3 does not change the third isospin component 
but it reverses the sign of f*^^) and t'-^\ Since C3——t'-"^\ 
we obtain from Eq. PT|) that 



-P- 



(48a) 
(48b) 



Consequently, in the absence of the explicit p-n mixing, 
the p-h density matrices have only the k=0 and 3 isospin 
components, while the p-p densities have only the k=l 
and 2 isospin components. 
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In the presence of an additional p-n exchange symme- 
try (charge-reversal transformation C of Eq. (1-5) multi- 
plied by i), defined as 

C2 = -m(f',07r0) = f(2), (49) 

only the fc=0 isospin component remains for p-h density 
matrix whereas the p-p has only a k=l no n- vanishing 
isospin component. In other words, in this case proton 
and neutron densities are equal to each other. 

B. Symmetries in the position-spin space 

In this section, we discuss transformation properties 

of the generalized density matrices TZ under D{a(3"f), P, 
and T. In the case of rotations, the general transfor- 
mation matrix ([28|) is given by d of Eq. (j43a|. Since 



the single-particle orbital angular- momentum operator is 
imaginary, l{r)= —l*{r), by applying ((3T|) to Eqs. (|43ap 
and (|44p. one obtains: 

d'^P-r^r's't', rst) = d""^^ {r' s' t' , rst), (50) 

i.e., the density matrices in both channels, p and p, trans- 
form under rotations in the same way, and the generalized 
rotation matrix has a simple form: 

W7)=(\ (51) 

When applying this symmetry operation to the gener- 
alized density matrix p4p . we need to use the Hcrmitian- 

conjugate matrix, (d°'^''\ : 



(d"'5'')^(r'sY,rsi) =<5(r'-r)5t-te-'^'-('^')e-''5'''(''')e-'"'-('-') x a+^{&,aPj) (52) 
and Pauli matrices that transform as vectors under the spin matrix a, i.e., 

d{&, af3j) (7° a+((T, af3j) = ^ aab{af3j) a'', (53) 

b 

for a,b — X, y, z, where the Cartesian rotation matrix reads 

( cos a cos /3 cos 7 — sin a sin 7, — cos a cos /3 sin 7 — sin a cos 7, cos a sin /3\ 
sin a cos/? cos 7 -I- cos a sin 7, — sin a cos /3 sin 7 -I- cos a cos 7, sin a sin /3 . (54) 
— sin [3 cos 7 , sin a sin fi , cos j3 j 

Similarly, the rotation matrix also rotates the position arguments, r and r', of the density matrices. Finally, we 
have 

p°{rst,r's't') = jpa{ar,sr')5ss'5tt' + \6ss'P{^r,sr') o fff/ + iso(ar, ar') • {acr)ss,5u' + \s{^r .s^r') ■ {a&)ss' oftt/,(55a) 
p°{rst,r's't') = \po{ar,ar')Sss'Stt' + jSss'P{3r,ar') o Tt^, + ^so{ar,ar') ■ {aajss'Stt' + js{ar,ar') ■ {a&)ss' oftt/.(55b) 

The inversion matrix (j43b|) is evidently real and symmetric, and it does not depend on a and r. Thus, we also 
have p{r' s't' , rst) — p{r' s' t' , r st) , and 

p^{rst,r's't') jPoi-r, -r')Sss'Stt' + jSss'P{-r, -r') o tu' + \sa{~r , -r') ■ &ss'Stt' + js{~r, -r') ■ ass' o ftt/, (56a) 
p^ {rst, r' s't') = \pa(-r, -r')6ss'5tt' + \5ss'P{-r, -r') o fu/-|-iso(-r, -r') • + is(-r, -r') • cTss' o ftt/. (56b) 

The time reversal T is an antilinear operation and has the form given in Eq. (1361) with the corresponding single- 
particle time-reversal matrix t given in Eq. I|43cp . It transforms all the position-dependent densities and the isospin 

Pauli matrices f to their complex conjugate partners and changes signs of cr. Therefore, t(r's't' , rst) — t{r's't' , rst), 
and the time-reversed densities are: 

p'^{rst,r's't') = lp*{r,r')Sss'Stt' + lSss'frir,r') of;,,-lsUr,r') ■ &ss'Stt'-lg*ir,r') ■ &SS' of:,,, (57a) 
^'^{rst,r's't') = \p*{r,r')6ss'Stt' + \Sss'P%r,r')of:,,-\s*{r,r')-ass'Su'-\i*{r,r')-ass'0^^^ (57b) 

I 



From Eqs. (j55bp . (j56bp . and (j57bp we conclude that the density matrix p in the p-p channel transforms under 
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all position-spin transformations considered here in the 
same way as the p-h density matrix p. In other words, 

the generalized transformation matrices of V and T have 

the same general structures as that of rotations, 2?, given 
in Eq. 

IV. SYMMETRY PROPERTIES OF DENSITIES 

We begin by recalling general symmetry properties of 
nonlocal densities, which have been given in I. Since the 
p-h density matrix (|8aP and the Pauli matrices are both 
Hermitian, all the p-h densities are Hermitian as well: 

Pk{r,r') = pUr',r), (58a) 
Sk{r,r') = sUr',r), (58b) 

for k — 0,1,2,3; hence, their real parts are symmetric, 
while the imaginary parts are antisymmetric, with re- 
spect to exchanging r and r'. 

Similarly, transformation properties of Pauli matrices 
under TC (|8bp imply that p-p densities are either sym- 
metric (scalar-isovector and vector-isoscalar) or antisym- 
metric (scalar-isoscalar and vector-isovector) under the 
transposition of their arguments, namely: 

Pk{r,r') = (-1)*'=+Vfe(r',r), (59a) 
h{r,r') = i-iy-hir^r), (59b) 

for fc = 0, 1, 2, 3, where = (isoscalars) and ti,2.3 — 1 
(isovectors). Equations ((SS]) and 159]) are fulfilled inde- 
pendently of any other symmetries conserved by the sys- 
tem; they result from definitions of density matrices p 
and p. 

A. Spherical symmetry 

1. Spherical and space-inversion symmetries 

Let us suppose that the generalized density 
matrix is invariant under the transformations of 
Eqs. (|43a|) and (|43b|l forming the full orthogonal group 
0(3)==P®SO(3)dSO(3), which is the direct product 
of the group of proper rotations SO (3) and the two- 
element group P of the space inversion. It means that 

iZ^ = TZ^ = HP^ = iZ, and these symmetries impose 
the following conditions on nonlocal densities: 



Pk{r,r') = pk{^zr,c,ar'), (60a) 

pk{r,r') = pfe(<^ar,<;ar'), (60b) 

Sk{r,r') = a+Sfc(<;ar,^ar'), (60c) 

Sk{r,r') = a+Sfe(^ar,^ar'), (60d) 



for k = 0,1,2,3, and for arbitrary Euler angles a, (3,^ 
that are arguments of the rotation matrix a (|44p . The 



factor is equal to -f 1 for rotations and —1 for improper 
rotations. 

The full 0(3) symmetry imposes quite strong condi- 
tions ([50]) on the nonlocal densities. Equations (|60ap 
and (j60bp tell us that, due to the Generalized Cay ley- 
Hamilton (GCH) theorem (see Appendix [X)) . scalar den- 
sities, pk and pk, depend on r and r' through rotational 
invariants r ■ r = r"^ , r' ■ r' = r'^, and r • r', i.e., 

Pk{r,r') = pk{r^,r ■r',r''^), (61a) 
Pkir.r') = pk{r\r-r\r'^), (61b) 

for k = 0,1,2,3. 

Similarly, from Eqs. (|60cp and (j60dp . we see that vec- 
tor densities, Sfc and Sk, are pseudovectors. At the same 
time, they are functions of two vectors, r and r' . The 
only pseudovector that can be constructed from two vec- 
tors is their vector product r x r' — therefore, all pseu- 
dovector densities, Sk and Sk, for fc = 0, 1, 2, 3, have the 
form 

Sk{r,r') i{r X r')sk{r'^ ,r ■ r' ,r''^), (62a) 
Sk{r,r') = (r X r')sfc(r^r • r',r'2). (62b) 

For the sake of convenience, in the definition (|62ap we 
have introduced the imaginary unit i. 

Due to the general symmetry properties ([55)1 and (|59p . 
scalar functions that define the nonlocal densities must 
obey the following conditions: 

Pk{r^,r -r'^r'^) = pl{r'^ ,r ■ r' ,r^), (63a) 

Pk{r\r-r\r'^) = {-lY''+'pk{r'^r ■ r' ,r^),{63h) 

and 

Sk{r'^,r -r^r'^) = sl{r'\r ■ r' ,r^), (64a) 

Sk{r^,r -r^r'^) = {-if^+^Skir'^ ,r ■ r' ,r^).{6Ah) 

This means that nonlocal p-h densities (|63aP and (|64ap 
are Hermitian, while nonlocal isoscalar and isovector p-p 
densities (j63bp and (j64bp are antisymmetric and sym- 
metric functions of r and r', respectively. 

The above symmetry properties of densities imply 
strong conditions on local densities. For instance, 
Eqs. (I6l])-(l64l) imply that 

Pfe(r) = pk{r,r) = pk{r) = pl{r), (65a) 

p{r)=p{r,r) = p{r), (65b) 

and 

Sk{r) = Sk{r,r) = 0, (66a) 

S(){r) = So{r,r) = 0, (66b) 

i.e., the scalar local p-h densities pk, and isovector scalar 
local p-p densities p, depend on the radial variable r only, 
Pk are real, and all vector local p-h densities Sk and the 
isoscalar vector local p-p density Sq vanish. At this point, 
we remind the reader that conditions (j59p imply that 
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local p-p densities po and s always vanish, irrespective 
of any other symmetries being imposed or not, see Table 
IV in I. 

All local derivative densities can be derived from 
Eqs. ((6T|) and ((62)) by using expressions for gradients of 
scalar functions, e.g.. 



dp 



dp 



Vp{r\vr'y^) = 2^r + ^^r', (67a) 
9P , + 2j4.r',(67b) 



V'p{r^,r ■ r\r'' 



d(r ■ r' 



9(r'2) 



which are linear combinations of vectors r and r' with 
scalar coefficients (again illustrating the GCH theorem) . 
In this way, all local derivative densities can be expressed 
through derivatives of the scalar functions pfc, pk, Sk, 
and Sk ■ Alternatively, one can employ the GCH theorem 
to build local scalar, pseudoscalar, vector, pseudovector, 
and symmetric-traceless-pseudotensor densities from the 
single position vector r, and we follow this path below. 

Pseudoscalar, pseudovector, and pseudotensor densi- 
ties cannot be built from the position vector r. There- 
fore, they must all vanish: 



Tk{r) 


= 0, 


(68a) 


To{r) 


- 0, 


(68b) 


Fkir) 


- 0, 


(68c) 


Mr) 


= 0, 


(68d) 


Jk{r) 


= 0, 


(68e) 


J{r) 


= 0, 


(68f) 


Ikabir) 


= 0, 


(68g) 


labir) 


= 0, 


(68h) 



for fc = 0, 1,2,3. 

The local scalar kinetic densities must have properties 
analogous to those of scalar densities in Eqs. ([SS]) . i.e., 

Tk{r) = Tk{r,r) = Tk{r) = T^{r), (69a) 
T{r)=f{r,r) = f{r), (69b) 

and all local vector densities must be proportional to r, 
i.e., 

jk{r) = jkr{r)er = jkrir)er, (70a) 
3o{r) = 3or{r)er, (70b) 
Jfc(r) = Jkr{r)er = Jl^{r)er, (70c) 

J(r) = Jr{r)er. (70d) 

where er- = - is the unit vector in radial direction. In 
addition, the radial components jkr and Jkr are real. 

Conditions on local densities, presented in this section, 
can be further restricted by imposing the time-reversal 
and/or p-n symmetries, see Table IV of I. The condi- 
tions on local time-even and proton-neutron-symmetric 
p-h and p-p densities are exactly those as in Ref. [33 |. 



whereupon properties of the p-p densities exactly mirror 
those of the p-h densities. Such a mirroring does not hold 
if time-reversal or p-n symmetries are broken. 

For a broken time-reversal symmetry and conserved 
p-n symmetry two modifications occur: (i) the isovec- 

tor p-p densities p(r), f{r), and Jr{r) become complex, 
while the p-h densities pk(r), Tk{r), and Jkrir) still re- 
main real [33|, and (ii) the current p-h density jkr{f) 
does not vanish. It is interesting to see that in this case 
the only non-zero time-odd density is the current density, 
i.e., spin polarizations are not allowed and only the flow 
of particles in the radial direction (a breathing mode) is 
permitted if the spherical and space- inversion symmetries 
are present. 

It is also interesting to see that the spherical and space- 
inversion symmetries impose very strong restrictions on 
the isoscalar pairing densities. Indeed, the isoscalar pair- 
ing density so{r) (j66bp must then vanish The only 
allowed isoscalar-pairing channel can be related to the 
p-p current density j'or (r) , which represents a radial flow 
of isoscalar pairs within a nucleus. Such a flow can be, in 
fact, non-zero either in the time-even or time-odd case, 
represented by ^{jor){r) or 3(jor)(?'), respectively. It 
corresponds to the situation in which p-n pairs locally 
change into the neutron-proton pairs, or vice versa, while 
the sum of densities thereof remains constant. 



2. Spherical symmetry alone 

Let us consider the unusual case of the SO (3) C 0(3) 
symmetry, in which the spherical symmetry of the gener- 
alized density matrix is conserved, whereas the space in- 
version symmetry is broken. As compared to the results 
presented in the previous section IIV A li here, proper- 
ties of scalar nonlocal densities, Eqs. (pT|) . remain the 
same. However, in Eqs. <r is always equal to -1-1, 
meaning that there is no difference between vectors and 
pseudovectors. Thus, vector nonlocal densities can now 
have structures that are richer than those of Eqs. (|62p . 
Indeed, due to the GCH theorem, these densities can be 
linear combinations of the pseudovector rxr' and vectors 
r and r'. 

Hence, vector densities in the p-h and p-p channels can 
now be presented in the form: 

Sfc(r,r') = i(r x r')sk{r^,r- r',r''^) (71a) 
+ rs'^{r^,r ■ r',r'^) + r's'k{r'^,r ■ r',r^), 

Sk{r,r') = (r X r')skir^ ,r ■ r' ,r'^) (71b) 
+ rs'kir^.r ■ r' ,r'^) + {~lY''r's'k{r'^ ,r ■ r' , r^), 

where scalar functions Sk and Sk obey previous conditions 
([64)) . while scalar functions s'f. and s'j. are arbitrary. 

Breaking the parity does not affect the local scalar den- 
sities in both channels and Eqs. ([55)1 and ([55)1 are still 
valid. The same is true for vector current densities of 
Eqs. (|70ap and (j70bp being gradients of scalar densities. 
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On the other hand, the local spin densities no longer van- 
ish. The isoscalar and isovector components of the p-h 
spin density are 



Sk{r) = Skr{r)er, 



(72) 



where their radial components s/cr are real for k = 
0,1,2,3. The p-p isoscalar spin density has a complex 
radial component and reads: 



So(r) = Sor{r)er. 



(73) 



Spin-kinetic and tensor-kinetic densities have the same 
structures as the spin densities of Eqs. ((72|) and ([73|) . 
namely 



Tfc(r)=T,!(r) = rfe,(r)e„ 
Fk{r) = F*^{r) = Ffc,(r)e, 



for all fc's and 



To(r) = To,(r)e, 



(74a) 
(74b) 



(75a) 
(75b) 



When the parity is not conserved, the tensor densities 
have, apart from the antisymmetric parts represented by 
vectors of Eqs. (|70cp and (|70d|) . also non- vanishing traces, 
namely. 



Jfc(r) = Jfe(r) = J,*(r), 
J(r) = J(r) 



(76a) 
(76b) 



and symmetric traceless parts have the following struc- 
ture: 

ifcab('') = 5iferr(^)Sa6 ~ ^iferr ('')Safc i ('''7a) 
lah{r) = \lkrr{r)Sat,, (77b) 

where the standard symmetric traceless tensor function 
of the space vector r is defined as 



— 3 "„ — 5ab- 



(78) 



3. Spherical symmetry - summary 



When the spherical 0(3) symmetry is conserved, all 
the local densities can be treated as fields depending on 
the 0(3) vector r (see Appendix E]). Then, vector r itself 
is the only one elementary vector and its length squared 

is the only one elementary scalar. The scalar densities 
are functions of r". The vector densities are pointing 
along r, and their radial components depend on only. 

The pseudoscalar, pseudovector, and pseudotensor 
densities all vanish. They can become non-vanishing 
when the parity is broken and only the rotational SO (3) 



symmetry remains conserved. Then, they have proper- 
ties of the scalar, vector, and tensor densities, respec- 
tively. The symmetric tensor densities are proportional 
to outer product r ® r. 

In general, the p-h densities are real whereas the p-p 
densities are complex. Additional symmetries (space in- 
version, time-reversal, p-n symmetry), when conserved, 
can cause additionally some isoscalars or isovector p-h 
densities to vanish, and those in the p-p channel become 
either purely real or purely imaginary, or vanish. In Ta- 
blesHlandim we list all the local S0(3)-invariant densities 
in cases when additional symmetries are broken or con- 
served. 



B. Axial symmetry 

While the spherical symmetry of the mean field is of- 
ten broken, the axial symmetry is usually conserved in 
the presence of time reversal. Here we discuss conse- 
quences of conserved symmetry of Eq. (|43ap with (3—0 
and 7=0, and the transformation of Eq. (|46cp . forming 
together group 0^^(2)=S^®SO-L(2)cO(3), which is the 
direct product of the rotations S0^(2) about the z-axis 
and the reflection Sz in the plane perpendicular to this 
axis. For rotations about the z-axis, the Cartesian rota- 
tion matrix (l54l) takes the form: 



cos a — sin a 
sin a cos a 
1 



a(aOO) 



a^(a) 
1 



(79) 



It is now convenient to decompose the position vectors 
as 



(80) 



where is a two-dimensional vector perpendicular to 
the z-axis and z is the z-component vector of r. 

Since -a^(a) = a-L(a + 7r), the 0^-^(2)= S^(8)SO^(2) 
symmetry of the generalized density matrix implies the 
following conditions for the nonlocal densities: 

Pk{r, r') = pfe(a^r_L + ^z, a-^r^ + qz'), (81a) 
Pk{r, r') = pfe(a^r_L + qz, a-^r^ + qz'), (81b) 
Sfe(r,r') = (r(a-L)+Sfci(a-Lr_L +<r2;,a^rl +<r2;') 
+SkAz^r^+qz,a^r'^ + qz'))e„ (81c) 
Sfe(r,r') = q{3.^)+Sk^{3.^r^ + qz,Si^r'^_ + qz') 
+Skz{z^r^+qz,s,^r'^ + qz')e^, (81d) 

where is the unit vector {z — ze^) and q is the sign 
of the determinant of the orthogonal 0^^(2) transforma- 
tion. 
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TABLE I: Properties of local particle-hole rotationally symmetric (SO(3)-invariant) densities, depending on the conserved (C) 
or broken (B) space-inversion (P), proton-neutron (p-n), or time-reversal (T) symmetries. Generic real, imaginary, or complex 
functions of the radial variable r are denoted by fnir), ,fi{r), or fcif), respectively. 
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Jl 
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n 


n 









n 
u 
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u 


J2 
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JR 
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,fij(r)er 








fR 
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TABLE II: Similar to Table |T] except for the particle-particle densities. 
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Axial and mirror symmetry 



In this case, Eqs. (|81ap and (|81bp imply that the scalar 
densities pk and pk depend on r and r' through the 
0^^(2) invariants z ■ z = z'^, z' ■ z' — z'^, z ■ z' — zz\ 
r±-r± — r^, r'^-r'j^ — r'^, and r±-r'j^. Apart from the in- 
variants, there are two 0^^(2) pseudoscalars: z-{r± xr^) 
and 2;'-(r^ xr^) = {zz' /z^)z-{r±xr'j_). The spin nonlo- 
cal densities have the following transformation properties 
under the 0^"'"(2). Their z-components are the 0^^{2) 
scalars. Their perpendicular components are the S0^(2) 
vectors and Sz pseudoscalars (i.e., they change sign under 
Sz)- Because the spin densities are the 0(3) pseudovec- 
tors, their components that are parallel to the z-axis 
should be linear combinations of r±xr'^, [z- {r± x r'^)]z, 
and [z' ■ {r± x r'j_)]z' , while the perpendicular compo- 
nents should be linear combinations of 2 x rj^, 2;' x r^, 
[z ■ {r± X r'j_)]r±, and [z' ■ {r± x r^)]r^. Consequently, 
the Hermitian nonlocal spin densities in the p-h channel 
should have the following structure: 

Skz{r,r')ez i{r^ x. r'jQkz(r,r') 

+ [z-{rx_ xr'^}]zQ'k,{r,r') 
^[z' ■[r^xr'^)]z'g'+{ry), (82a) 
SkAr.,r') = {z X r±)gk±{r,r') 

+ iz'xr'^)gtAr,r') 

+ [z ■ (r_L X r'j^)]r±g'^^^{r,r') 

-[^'.(r^xr^)Ke;\(r,r'),(82b) 

where gkz, dkz^ 8k±, and p'^,^ are scalar functions. The 
pseudovector nonlocal densities in the p-p channel are ei- 
ther symmetric or antisymmetric in r and r' . Therefore, 
we have 

Soz{r,r')ez = {r± x r'j_)goz{r,r') 

+ [z ■ {r± X r'^)]zg'oz{r,r') 
-[z' ■{r^xr'^)]z'g'oz{r\r), (83a) 
Sz{r,r')ez = {r± x r'j_)gz{r,r') 

+ [z ■ (r^ X r'^)]z^'^{r,r') 
+ [z' ■ir^xr'^)]z'l{r',r), (83b) 
so±{r,r') = (z X r±)go±{r,r') 

+ {z' X r^)^o±(r',r) 
+ [z ■ (r_L X r'^)]r^g'„j_{r,r') 
-[z' ■{r^xr'^)y^g'oAr\r),i83c) 
s±{r,r') = (z X r±)g±{r,r') 

-{z' X r'A^Ar',r) 

+ [z ■ (r_L X r'j^)]r^g'j_{r, r') 

+ [z' ■ {r^ X r'A]r'J'Ar',r), (83d) 

where ^oJ_ , q'oi_ i Qi- ^ q'i_ 1 q'qz 1 and are arbitrary scalar 
functions, while g^z and gz are antisymmetric and sym- 
metric scalars, respectively. 



The local scalar p-h and p-p densities are functions of 
two invariants, and r\. The particle density pkiz'^, t'^) 
for k = 0,1,2,3 is real whereas the pairing density 
p{z'^,r']_) is complex in general; hence, 

Pk{r) = pk{r,r) = pk{z\rl) ^ pl{z^ ,rl), (84a) 
^{r)=^ {r,r) = ^{z\rl). (84b) 

Differential operators Vz and V_l have the same trans- 
formation properties under the 0^"'"(2) transformations 
as z and r±, respectively. That is, Vz and z are Sz 
pseudo-invariants and SO-*" (2) invariants, whereas Vj_ 
and r± are Sz invariants and SO-"" (2) vectors. Coun- 
terparts of Eqs. (I67ap and (I67bp for gradients of scalar 
functions are now linear combinations of vectors z, z' , 
r±, and r'j_ with scalar coefficients: 

(Vz + Vi)p(z2, zz', z'2, rl,r^- r'^,r'l) 



= 2^z^ 
9(z2) 

„ dp 



dp ^, 



d{zz' 



dp 



' d{r'\) ' ^ ' 9(rj_ • t^) 
(V; + Vl)p(z2, zz', z'\ rl,r^- r^, r'l 



TT^l- (85a) 



-I- 



d{zz') a(z'2)' 

dp dp 



d{r± ■ r'j_) d{r'\ 
Both terms of the operator: 

V ■ V = Vz ■ vt + V I ■ V'l 



■r^L- (85b) 



(86) 



are 0^-'"(2) scalars. Therefore, the local scalar kinetic 
densities: 

Tk{z',rl) = [(Vz- V'z + Vi- Vl)pfc(r,r')].=w(87a) 
f(z2,ri) = [(Vz- V; + Vi- Vl)/5(r,r')].=w,(87b) 

can be expressed as sums of the two O^-*" (2) scalars. Since 
the operator (j86p is Hermitian, again the p-h densities are 
real and the p-p ones are complex. 

It is seen from Eqs. (|82|) and (|83l that the spin densities 
in both channels are parallel to the vector product zxr± 
and take the form: 



Ba) 



Sk{r) = gk±{z'^,r1^){z X r±), 
So = gQ±{z'^,r]_){z x r±), 

with real gki_ and complex go±. Applying Eq. (l85|) to 
Eqs. ([5^ and (|83p . we find that the spin-kinetic and 
tensor-kinetic densities in both channels can be written 
as: 

Tfc(r) = [(Vz- V; + Vi- Vl)sfc(r,r')]r=r' 

= Mz^rDizxr^, (89a) 

To(r) = [(Vz- V'z + Vi- Vl)s^o(r,r')].=.' 

= Mz^,rl)izxrA, (89b) 
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and 

Fkir) = ^[(V, + V^)((V', + Vl).s,.(r,0) 

+ (V; + Vl)((V, + Vi) • Sk{r, r'))]r=r' 
= ^k{z^,rl){z X r^), (90a) 

Foir) - i[(V, + V^)((V; + Vl).so(r,r')) 

+ (V; + V1)((V, + Vi) • so(r,r'))].=w 
= ^o(2^ri)(^ xr^). (90b) 

The scalar functions and are real whereas i?o and 
(/?o are, in general, complex. Equations (I85p imply that 
the vector current densities in both channels are spanned 
by vectors z and r±, and read 

jfel^") = ikz{z^,r'^±)z + Lk±{z'^,r'i)r±, (91a) 
Jo(t-) = ioz{z^,ri)z + Lo±{z^,rl)r±, (91b) 

where t^z, tfcx, lqz, and to± are scalar functions. The p-h 
currents are real and the p-p current is complex. 

The spin-current pseudotensor densities can be decom- 
posed according to: 

J,(r) = J^.(r)+J^(r) + Jf (r), (92a) 

J(r) = y{r) + y{r) + y^{r), (92b) 

where 



Jfe(r) 




- V'J Sfez(r,r')ez]r=r 


, (93a) 


J^(r) 




- Vl) «> Sfei(r,r')]r=^/ 


(93b) 


(r) 




- V'J®Sfei(r,r') 






+ (V± - 


Vl) ® Sfcz(r,r')ez]^=w 


, (93c) 


J^(r) 


= 


- V'J®S4r,r')ez]^=w, 


(93d) 


J^(r) 




- V'L)®Si(r,r')]r=r', 


(93e) 


P^(r) 




- V^) S^(r,r') 






+(Vi- 


Vl) ® Sz(r,r')e2]^=^/. 


(93f) 



The traces of the spin-current pseudotensors: 

Mr) = Jl{r) + Jt{r), (94a) 
J(r) = r{r) + J^{r), (94b) 

where 









^kziy ^'^ )^z\r—r' 


(95a) 


Jt{r) = 




-VI) 




(95b) 


J\r) = 








(95c) 


J^(r) = 




-VI) 


■ s±{r,r')]r=r', 


(95d) 



are the scalar products of gradient operators and spin 

densities. Since Jk and J are sums of 0^^(2) pseu- 
doscalars, they cannot be constructed from two vectors 
z and r±; hence, all spin currents (|95p must vanish. On 
the other hand, the 0(3) vectors coming from the anti- 
symmetric parts of spin-current pseudotensors: 

Mr) = ^[(V,- V'Jx8,^(r,r') 

+ (Vi - VY) X Sk{r,r')]r^r', (96a) 

Jir) = :^[(V,-V;)xl^(r,r') 

+ (Vi-Vl)xl(r,r')].=w, (96b) 

do not vanish. They can be decomposed in the same way 
as the current vectors (|9ip : 

Jk{r) = Vkz{z'^,r'i)z + Vk±{z'^,r'i)r±, (97a) 
Jir) = Sziz^rl)z + 6^{z^,rl)r^. (97b) 

The scalar functions Vkz and Vk± are real whereas ijz and 
v± are complex. Finally, the traceless symmetric parts 
of the spin-current densities (|92p are: 

llir) = 0, (98a) 

r(r) = 0, (98b) 

J;^(r) = Ki{z',rl)ir^^iz^r^), (98c) 

//(^) = 4^(z^ri)( ^0(^xr^) ), (98d) 

V{r) = ^^(z^ri)( r^ (z x r_L) ), (98e) 

J"^(r) = ^"^(^',ri)( z0(2 xr^) ), (98f) 

where 

r_L (g) (2; X r^) = ^[r± ® (z X r±) 

+ (2; X r^) (g)r^] , (99a) 
z (g) (2; X r±) = i [z (g) (2; X r^) 

+ (2xri)g)2] (99b) 

are the symmetrized outer products of the vector product 
zxr±, and vectors r± and z, respectively. As usual, the 
scalar functions Kk are real, whereas k are complex. 

2. Axial symmetry alone 

From the above discussion we see that the mirror sym- 
metry imposes quite strong conditions on the nonlocal 
and local density functions. When the generalized den- 
sity matrix is invariant only under the SO^(2)cSO(3) 
group of rotations about the z-axis, the transformation 
rules of Eqs. ([5T|) with <r = -f 1 are fulfilled. This means 
that coordinate = z is the S0^(2) invariant, and there 
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is no difference between the SO (3) pseudovectors and vec- 
tors. The scalar nonlocal densities are thus functions of 
rj_, 2', r^, and r \_ ■ r'^. 

The spin densities are S0-'-(2) vectors and can take 
more general forms than those of Eqs. ([82|) and ([83|) : 
namely, the Hermitian spin densities in the p-h channel 
are: 

Skz{r,r')ez = i{r± x r'^)gkz{r,r') 

+zgUr,r') + z'g'C,{r\r), (100a) 
Sfc_L(r,r') = {zx ri_)gk±{r,r') + {z' x r'^)gl^{r' ,r) 
+r^g'kAr, r') + r'^g'CAr' , r), (100b) 

where gkz is Hermitian and g'j,^, gk±, and g'j,^ are arbi- 
trary scalars. 

The spin nonlocal densities in the p-p channel are ei- 
ther symmetric or antisymmetric in r and r': 

Soz(r,r')e^ = {r± x r^)^o^(r,r') 

+zg'oAr,r') + z'g',,ir',r), (101a) 
s^(r,r')e^_ = (r_L x r^)^^(r,r') 

+2;ei(r,r')-^'g*Ur-',r), (101b) 
So±(r,r') = (2; X ri_)gQ^{r,r') + (2;' x r^)^o±(T'', r), 

+r^g'^^{ry) + r'^g'^^{r',r), (101c) 
s_L(r,r') = (2; X r^)g_L(r',T') - {z' x r^)p^(r',r) 

-Krill(r,r')-rYll(r',r), (lOld) 

where g^z is antisymmetric, gz is symmetric, and gg^, 
Po-L, ^o_lj Pzj Pij ^'^d g'^ are arbitrary complex scalar 
functions. 

In the case of broken mirror symmetry, the real local 
p-h densities Pfc(2, 7'j_) = p%{z^ rj^) and the complex local 
isovector p-p density /9(z,rj_) are functions of the two 
cylindrical coordinates, z and r±, i.e., they depend on the 
sign of z. The same is true for the scalar kinetic densities 
([57)) . i.e., the real p-h densities Tk{z, r_i) = t^{z, r±) and 
complex p-p isovector kinetic density -r(z,r^). 

At this point, it becomes convenient to use the cylindri- 
cal coordinates r±, 0, and z, and the corresponding unit 
vectors e±, e^, and Cz. The p-h isoscalar and isovec- 
tor real spin densities and the p-p isoscalar complex spin 
density are: 

Sfc(r) = Skr^{z,rj_)ej_ 

+ Sk,i,{z,rj_)e^ + Skziz,rj_)ez (102) 

for fc = 0, 1, 2, 3, and 

so{r) = SQr^{z,r±)e± 

+ SQ^{z,r±)e^ + Soz{z,r±)ez. (103) 

All differential local densities can be calculated by using 



the gradient formulae: 



(Vz + Vj^)p(z, z', r^,r^- r'^,r\ 
dp dp 



dz ^ dr± d{r± ■ r'^) 
{V', + V'^)p{z,z\r^,r^-r'^,r'l) 

oz' o(r± ■ r^) ar^ 

The spin-kinetic densities ([55)1 have all non-vanishing 
cylindrical components: 

Tfe(r) = Tkr^{z,r^)ei_ 

+ Tk^{z,r±)e^ + Tkz{z,r±)ez, (105a) 

fo(r) = fo,^(z,ri)ei 

+ fo4,{z,ri_)e^ +fa,{z,r^_)ez. (105b) 

As usual, the p-h spin-kinetic densities are real and the 
p-p ones are complex. 

The tensor-kinetic densities (|90ap and (j90bp have the 
same structure as the vectors (|105p : 

Fk{r) = Fkr^iz,r±_)e± 

+ Fk^{z,r±)e^ + Fkz{z,r±)ez, (106a) 

Fo{r) = FQr^{z,r±_)e±_ 

+ Fo4,iz,r±_)e^ + Faz{z,ri_)ez. (106b) 

The current densities, being proportional to the gradients 
of scalar functions, have only the rj_- and z-components: 

jk{r) ^ jkz{z,r±)e:,+jkr^{z,r±)e±, (107a) 
jo{r) ^ joz{z,r±)e:,+jor^{z,r±)e±. (107b) 

Finally, due to the breaking of mirror symmetry, the 
spin-current densities can have rich structures. The 
traces of spin-current tensors ([55]) do not vanish, and 
they decompose into the sums of two S0-'-(2) scalars: 

Jk{r) = Jl{z,r^) + J^{z,r^), (108a) 

J(r) = J\z,r^) + J^{z,rj_). (108b) 



The antisymmetric parts (|96p of the spin-current tensors 
(1931) form the vectors with nonzero transverse compo- 
nents: 

Jfe(r) = J^(r) = Jkr^{z,r±)e± 
+Jk4>iz,rj_)e^ + Jkziz,r±)ez, (109a) 

J(r) = Jr^{z,r±)e± 

+J^{z,r±_)e^ + Jz{z,r±)ez. (109b) 
The symmetric traceless parts of the tensors (|M)) vanish: 



-kzz 



= 0, 
= 0. 



(UOa) 
(110b) 
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The remaining traceless symmetric tensors are: 

—kab :^kah 

= Kfe_L±(z,r±)Pab + Jfc±_L(2:''-L)Sa6, (Ula) 

L = i^A^,r^)Pi,+l^^{z,r^)Si, (111b) 
for a,b — X, y, where 

T I „T I J, 

(112) 



9.ab — ^ 2 "if' 



is the standard symmetric traceless S0^(2) tensor and 

^ab = ol % ir±a£bzc + r±b£azc)r±c (113) 



2|z|ri 



is a symmetric pseudotensor. The tensor is trace- 
less by definition. Its symmetric part has the following 
structure: 



iLz = Ikaz = \<LkzAz,r±)Pl, +Jl,^^{z,ri_)Sl^ (114) 

for a — x,y, where non-vanishing components of the 
normalized and symmetrized outer products z (E>r± and 
z (g) (z X r±) are: 



and 



Similarly, 



r± \z\ 



iaz = R2i(^:,f±)Pa^ +i±(z,r_L)S^^ 



(115) 



(116) 



(117) 



All components of the spin-current tensor in the p-h 
channel are real whereas those in the p-p channel are 
complex. 



3. Axial symmetry - summary 

In the case of the axial 0^^(2) symmetry, the position 
vector r can be decomposed ([50)1 into two vectors z and 
r_L having different transformation properties under ro- 
tations and mirror rotations about the z-axis. Vector 
is the SO"*- (2) vector, whereas z is the 0^^(2) pseudo- 
invariant. There are two 0^^{2) scalars: and r\, and 
all the local scalar densities are functions thereof. In this 
study, we are not concerned with the question whether 
the densities are analytical functions of the invariants. 
Therefore, it does not matter whether the argument of 
densities is r\ or just r±_. 

The local vector densities are linear combinations of 
vectors z and r±_ with scalar coefficients. The pseu- 
dovector densities are proportional to vector product 
z X r± and thus have the azimuthal direction. While 



it is not possible to build pseudoscalar densities from ele- 
mentary vectors z and r± , pseudotensor densities can be 
constructed. The symmetric pseudotensor densities are 
linear combinations of the symmetrized outer products 
r±® [z X r±) and z ® {z x r^). The symmetry proper- 
ties of local axially symmetric (0^^(2)-invariant) densi- 
ties are listed in Tables IIIII and IIVI 

If the mirror symmetry is broken, the two S0'''(2) in- 
variants are z and r±. The scalar and pseudoscalar den- 
sities are now functions thereof. The vector and pseu- 
dovector densities now have non-vanishing components 
along all the three vectors, ej^, e^, and e^- The trace- 
less symmetric tensor densities are linear combinations 
of the pseudotensors and P^"*", and tensors and 
S^^. Properties of the S0^(2)-invariant local densities 
are listed in Tables |V] and |Vll 



C. Symmetry D2h 

The identity, inversion, three signatures, three sim- 
plexes, and their negative partners form the symmetry 
group D^jj . The Cartesian rotation matrices of Eq. ((54|) 
for the identity (labeled by u) and signature operations 
are all diagonal: 



a° 

a;; 



for a = u, X, y, z. They can be written explicitly as 



(118) 



a" = a(OOO) 
a^ = a(07r7r) 
= a(07rO) 
a^ = a(7rOO) 



1 




1 


)■ 


1^ 




1 





-1 








-1 


-1 





1 





^ 


-1 


-1 





-1 








1 



(119a) 
(119b) 
(119c) 
(119d) 



Let us suppose that the generalized density matrix is in- 
variant under the generalized transformation matrix U"", 
(a = u, X, y, z), belonging to the D^j^ of transformations. 
According to Eqs. (|60p . the transformation rules for non- 
local densities are: 



Pkir,r') 
Pk{r,r') 
Skb{r,r') 
Skb{r,r') 



Pfe(^aV,^aV), 
Pfe(^aV,^aV), 
a^Sfeb(^aV,^aV') 
a^Sfeb(caV,<;aV') 



(120a) 
(120b) 
(120c) 
(120d) 
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TABLE III: Properties of the local axial and mirror-symmetric (0^"'"(2)-invariant) particle- hole densities, depending on the 
conserved (C) or broken (B) proton-neutron (p-n), and time-reversal (T) symmetries. The z-simplex (Sz) symmetry is conserved. 
Vector and pseudovector densities can be expanded in a basis of three vectors: r±, z and z x r±. The symmetrized outer 
products r± iSi {z X r±) and z {z x rj_) form a basis for the pseudotensor densities. The expansion coefficients are real, 

2 



imaginary, and complex functions of the two (2) scalars: z 
denoted by /^(z^, r^^), //(2'^,r'^), and fc{z^,r^), respectively. 



and r^. Generic real, imaginary, and complex coefficients are 



TT 

B 
B 



conserved (C) or broken (B) 



TT 

C 
C 



basis 



coetticients 



p0,3 


I 










fc 

JR 


\Z 5^"xj 


JR\Z , 1 \) 
JR\Z ,r^) 


J 




fa(z'^ \ 
J R\^ j'±) 


Pl 


1 










jR 


n 

u 




n 

u 


p2 


\ 










JR 




Q 



u 




Q 


"^"0.3 


\ 










JR 


/ 2 2 \ 

)^xj 


J R{z , r ±) 






J R\Z , 1 ±) 


'^1 


\ 










JR 


',2 2 N 


JR(Z ,r^j 








u 


T2 


1 
1 










JR 


n 
u 


n 
u 




n 
u 


Jo,3 





























Ji 





























J2 





























So, 3 


z 


X r± 








JR 


r 2 2 \ 


n 

u 


JR\Z 




n 

u 


Si 


z 


X r±_ 








fR 














S2 


z 


X r± 








fR 




)R{z\rl) 










To. 3 


z 


X r± 








fR 


(2 2 \ 

z ,r , ) 
r 2 2 \ 

(2 2 \ 

yZ ,r A 







rl) 





Ti 


z 


X r± 








fR 












T2 


z 


X r± 








fR 


fR{z\rl) 












z 


X r± 








fR 


1 2 2 \ 







rl) 





Fi 


z 


X r± 








fR 


'z^ ) 












F2 


z 


X r± 








fR 


fR{z\rl) 










jo, 3 


r 


L,Z 








fR 


[z^rl) 





fR{z^ 


rl) 





jl 


r_ 


L,Z 








fR 


1 2 2 \ 
z ,r±) 













32 


r_ 


L,Z 








fR 


1 2 2 \ 
z ,r±) 

[z\rl) 


fRizyi) 









Jo, 3 


r 


L,Z 








fR 


fR{z\ri) 


/«(^' 


rl) 


fR(zlrl) 


Jl 


r 


L,Z 








fR 


[z\rl) 


fR{z',rl) 










J2 


r_ 


L,Z 








fR 


1 2 2 \ 
z ,r±) 











Jo. 3 


r_ 


L®{Z 


X rj 


-),Z(i 


5 (2: X rx) 


fR 


r 2 2 \ 

^z%r\_) 


fR{z\rl) 




rl) 


fR(z\rl) 


ll 


r 


L®{Z 


X rj 


J, 2« 


5 (2: X rx) 


fR 


(2 2 \ 

^z%r\_) 


fR{z\rl) 










I2 


r 


L®{Z 


X rj 


^,z(^ 


5 (z X rx) 


fR 


(2 2 \ 

^z ,rY) 






































for any a — u,x,y, z and b = x,y, z; <; = +1 for identity 
and signature operations and = — 1 for inversion and 
simplex operations. 

The transformations (|120p constitute symmetry condi- 
tions for densities under changes of signs of their argu- 
ments; hence, they relate the values of densities between 
different regions of space. For the p-h densities, this prob- 
lem has been extensively discussed in Refs. [22|, [231 • In 
the following, we extend the previous discussion to the 
pairing channel and provide general expressions for the 
p-p densities. 

The p-h local density matrices fulfill the following sym- 
metry conditions: 

Pk{r) = pKr)=pfe(?aV), (121a) 
Skb{r) = sl,{r) = a?Sfc,(^aV) (121b) 

for A; = 0, 1, 2, 3, a = u, x, y, z and b — x,y, z. The anal- 
ogous expression for the p-p densities are: 



p(r) = p(<;aV), 
So6(r) = a;^Sob(<;a''r). 



(122a) 
(122b) 



The differential local densities in both channels should 
now be classified according to irreducible representations 
of the point group D2h [l^l- Therefore, the vector no- 
tation used so far is no longer useful. Instead, in this 
section we rely on definitions p^ . where the Cartesian 
components are explicitly shown. 

The p-h isoscalar and isovector kinetic densities (|17ap 
transform according to: 



Tfebc(r-) = T4,(r) = a;ja^Tfe;„(^aV) 



(123) 



The complex isovector component of p-p kinetic densities 
(I17b[) obeys analogous symmetry conditions: 



hc{r) = a?a^Tf,,(^aV) 



(124) 



The symmetries of the kinetic scalar densities ([TT|) can 
be obtained from Eqs. (|18ap and (jl8b[) . 

The spin-kinetic and tensor-kinetic densities (|17cp and 
(I17dp fulfill the symmetry conditions: 

Tu^>rAr) = T4,,(r) = a«a«a«Tfc,,d(?aV), (125a) 



To6crf(r-) = a^a:;aJ^rob,d('?aV) 



(125b) 



17 



TABLE IV: Similar to Table Hill except for the particle-particle densities. 



symmetry 








conserved (C) 


or broken (B) 










C 


C 


C 


c 


p-n 






B 


B 


C 


c 


r 






B 


C 


B 


c 




basis 






coefficients 




pi 

P2 
P3 
Tl 
T2 
T3 


1 

1 
i 

-1 
i 

1 
1 

1 




/c(^^rl) 
/c(2 ,rx) 

/. /2 2 \ 

jc(z ,rj_) 
jc(z ,rxj 

./Cl,^ , ' x^ 

fc(z^ ) 


fi{z\rl) 
fR{z\rl) 

jRyz , 1 A_) 

fjiz^ r^) 
J , / _l; 


/c(^^rl) 
/c(^^ri) 
U 

JC\Z , I j_) 




fR{z\ri) 

f (',2 „2 \ 
U 

J ' x,i 




Jl 



















J2 



















J3 


















So 


z X r± 




fciz^ri) 


fi{z\rl) 








To 


z X r± 




fc{z',rl) 


fiiz^rl) 








Fo 


z X r± 




fc{z\rl) 


h{z\rl) 








30 
Jl 


r±,z 
r±,z 




fc{z\rl) 
jc(z ,rj_) 


fi{z\rl) 
fR(z\rl) 
fi{z\rl) 




/c(^^ri) 




Jr(z ,r±) 
Mz',rl) 


J2 


r±,z 




fc{z\rl) 


/c(^^ri) 


J3 


r±,z 




fc{z\rl) 


fR{z',rl) 








ii 


(g) (z X rj_), z (J 


5 (z X rx) 


fc{z\rl) 


fR{z\rl) 


/c(^^ri) 


fR{z\ rl) 


I2 


rx (g) (z X r±), z $ 


5 (z X rx) 


fc(z\rl) 


fiiz',rl) 


/c(.^ri) 


Mz^rl) 


I3 


rx ® (z X rx), z $ 


5 (z X rx) 


fc{z\rl) 


fRiz\rl) 























The transformation properties of the spin-kinetic densi- 
ties and tensor-kinetic densities (fT5|) are obtained by 
contractions defined by Eqs. (|18c|) - (|18fp . 

Finally, the symmetry conditions of the current and 
the spin-current densities in both channels are: 

jkbir) = jl,(r) = ^a^j,b(^aV), (126a) 
jobir) - ?a^job(?aV, (126b) 
hbcir) = JI.,,(r)=^a^a^Jfc,,(^aV), (126c) 

L(r) = ^a^a^L(';aV). (126d) 



What are the differences in the symmetry properties of 
the transition density matrices from those of the density 
matrices discussed above? First of all, the p-h transition 
density matrix is not Hermitian: 



{rst,r's't') = («'i|a+^,j,a^^t|*2) 

7^ (*2K,,,t,ar.t|*i). (129) 



P 



Consequently, Eqs. ([58]) for the p-h nonlocal transition 
densities are not fulfilled. On the other hand, the an- 
tisymmetry property of the p-p density matrix is pre- 
served: 



V. SYMMETRIES OF MULTI-REFERENCE 
TRANSITION DENSITIES 

In analogy to Eqs. ^ and the transition p-h and 
p-p density matrices are defined, respectively, as 



s(t) 



{rst.r's't') = (^'2|a+^,t,a^^t|*i), 



(127) 



As't'{^2\ar'-s'-t'arst\'fl),{l28) 



where and |5'2) are two different independent- 

quasiparticle states. The corresponding spin-isospin 
scalar and vector transition densities p]^\r,r') and 



s^^\r, r'), and p)^'{r, r') and s'"^'{r, r') with /c = 0, . . . ,3 
are defined by the relations analogous to Eqs. ([9]). The 
local transition densities are defined in an identical way 
as the local densities Eqs. PIH) - (IT5|) and are denoted with 
the same respective symbols but with the superscript {t) . 



(*)/ 



=(t) 



{rst,r's't') = -16ss'tt'/B(*)(r' - s' -t',r- 



s - t). 
(130) 

Hence, the p-p nonlocal transition densities obey rela- 
tions ((59)) . 

The transformation rules for the transition matrices 
under the single-particle unitary (and antiunitary) trans- 
formations U [Uk) follow the transformation rules P7)) 
for the creation and annihilation operators and are given 
by Eqs. (|5D|). namely: 



't)U _ 
(t)U 



(131a) 
(131b) 



Therefore, the discussion of density matrix symmetries 
presented in Sec. IIVI applies to transition densities with 
the only difference being that the p-h transition densi- 
ties are, in general, complex, unless the time-reversal in- 
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TABLE V: Properties of local axially symmetric (S0^(2)-invariant) particle- hole densities, depending on the conserved (C) 
or broken (B) p-n and time-reversal symmetries. The z-simplex (Sz) symmetry is broken. The vector, pseudovector, or 
pseudotensor densities can be expanded in the vector (ex, e^, Ez) or tensor S^, S^^) bases. The expansion coefficients 

are real, imaginary, and complex functions of cylindrical coordinates rx and z. Generic real, imaginary, and complex coefficients 
are denoted by fR{z,rj_), fi{z,r±), or fc{z,r±), respectively. 



symmetry 












conserved (C) 


or broken (B) 












B 






B 






B 




B 




p-n 






B 






B 






C 




C 




T 






B 






C 






B 




C 






basis 














coefficients 








po,a 


1 




/fl 


z 




fR 




r 1 




; 

rx) 


fR{z 


r±) 


Pi 


1 




}r 


z 




fR 




T 1 1 












p2 


1 




Jr 


z 


rx' 



















To, 3 


1 




Ir 


z 


rx. 


fR 




r 1 ) 


fnlz 
.1 


rx ) 


fR^z 


r±) 


Tl 


1 




fR 


z 




fR 




T 1 1 












T2 


1 




fR 


z 


rx' 



















Jo,3 


1 




fR 


z 


r± 


fR 






fR{z 


rx J 


fRiz 


r±) 


Jl 


1 




fR 


z 


^x, 


fR 




rxj 












J2 


1 




fR 


z 





















So, 3 


ex, 60, 


6z 


fR 


^z 











fR{z 


rx) 







Si 


ex, e^. 


6z 


fR 


' z 





















S2 


ex, e^. 


6z 


fR 


z 




Jr{z 


rx) 












To, 3 


ex, 60, 


6z 


fR 


z 


r±] 









fR{z 


rx) 







-t 1 


ex, 60, 


6z 


J- 

fR 


z 


r±. 














u 




T2 


ex, 60, 


62 


fR 


z 


rx. 




rx) 












^0,3 


ex, 60, 


6z 


fR 


z 


^X; 









fR{z 


rx) 







Fi 


ex, 60, 


6z 


fR 


z 


rx' 



















F2 


ex, 60, 


6z 


fR 


z 


rx. 


fR{z 


rx) 












jo, 3 


ex, 60, 


6z 


fR 


z 


^X; 









fR(z 


rx) 







jl 


ex, 60, 


6z 


fR 


z 





















32 


ex, 60, 


6z 


fR 


z 




fR 




rx) 












J0.3 


ex, 60, 


6z 


fR 


z 




fR 


i 


rx) 


fR{z 


rx) 


fRiz 


r±) 


Jl 


ex, 60, 


6z 


fR 


z 




fR 




rx) 












J2 


ex, 60, 


62 


fR 


z 


r±] 



















=io,3 


pX p. 




fR 


z 


r±] 




rx) 


fR{z 


rx) 


fRiz 




il 


pX p. 




fR 




r±] 


fR{z 


rx) 












J2 


pX pz 




fR 


> 


r±] 




















variance introduces some restrictions. The time-reversal- 
invariant non-local p-h transition densities obey the fol- 
lowing relations: 

pi'\r.r')=p^^>{r,r'), (132a) 
4*^(r,r') = -4*^*(r,r') (132b) 
for fc = 0, 1, 3, 



and 



p^^\r,r') = -p'^>*{r,r'), (132c) 
4*)(r,r') = 4*^*(r,0, (132d) 

which means that some transition densities are real while 
others are purely imaginary. Symmetry properties of the 
local p-h transition densities are catalogued in Tables 
IVIIl IVIlTl and IIXI For the local p-p transition densities, 
the results of Tables IH |TVl and |Vl] apply. 



VI. DISCUSSION OF PAIRING CHANNELS 
AND EXAMPLES OF PREVIOUS APPROACHES 

The form of the most general EDF that is quadratic in 
local isoscalar and isovector densities has been proposed 



in Ref. I, where the expressions for the p-h and p-p mean 
fields can be found. In current applications, pairing in- 
teraction is often approximated by the zero-range pairing 

force HIlEIsi,!!!, 



^pair(r-,r') = /, 



pair 



{r)S{r — r') 



(133) 



where the density-dependent form factor reads 

>o(r)" 



/pair(r)-l^O 1+xoP' 



Pc 



(134) 

and J"^ is the usual spin-exchange operator. When only 
the isovector pairing is studied, the exchange parame- 
ters a:o and X3 are usually set to zero. However, in the 
general case of coexisting isoscalar and isovector pairing 
correlations, nonzero values of xq and have to be used. 

In Ref. |38j . the density-independent, zero-range pair- 
ing force, 

Vp.ir{r,r') = [Po^i'' - r')+p^k'Sir - r')k] Hts, 



TS 



(135) 

has been employed to study the interplay between 
isoscalar and isovector pairing within an axially symmet- 
ric HF+BCS scheme. In Eq. tlrs stands for the 
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TABLE VI; Similar to Table FVl except for the particle-particle densities. 



symmetry 








conserved (C) 


or broken (B) 










B 


B 


B 


B 


p-n 






B 


B 


C 


C 


r 






B 


C 


B 


C 




basis 




coefficients 




pi 


1 




fc{z,r±) 


fR{z,r±) 


fciz,r_i_) 


fR{z,r±) 


P2 


1 




fc{z,r±) 


fi{z,r±) 


fc{z,r_L) 


fi{z,r_L) 


P3 


1 




fc{z,r±) 










fl 


1 
1 




.fc{z,r±) 


fniz, r±) 


fc{z:, rj_) 


/i?(z,rx) 


T2 


\ 




fc{z, r±) 


fi{z, r±) 


Ti^ ( 'y or" 1 \ 

JC\Z^ f _L) 


fi{z, r±) 


T3 


1 




fc{z, r±) 


fR{z^ r±) 








Jl 


1 




fn ( z V \ \ 


fftiz.r 1 ) 


fr( z.r \ ) 

J ^ 7 ' -L J 


fr? ( Z T ] ^ 

J R\^-> ' ^ ) 


J2 


1 




fn ( z V \ \ 


fi{z,r_L) 


fc{z,rj_) 


friz T ] \ 


,h 


1 




fn ( z V ] ) 


fR{z,r±) 








So 


e_i 




fc{z, r±) 


.fi{z,rj_) 








To 


e_i 


^<f>1 


fc{z, r±) 


fi{z,r±) 










e_i 


_ 7 6(/> ; ^2 


fc{z, r±) 


fiiz,r±) 








30 


e_i 


_ 7 ^(/>! 


fc{z,rj_) 


fi{z,r±) 








Jl 


e_i 


_ 7 ^1^5 


fc{z,rj_) 


fR(z,r±) 


fciz,ri_) 


fR{z,rj_) 


J2 


e_i 


_ J ^(/!>5 ^2 


fc{z,r±) 


fiiz,r±) 


fc{z,ri_) 


fiiz,r±) 


J3 


e_i 


-■) ^<j)t ^2 


fc{z,r±) 


fR(z,r±) 








ii 






fc{z,r^) 


fR{z,rj_) 


fc{z,rj_) 


fR{z,r±) 


I2 






fc{z,r^) 


fiiz,r±) 


fc{z,rj_) 


fiiz,ri_) 


I3 


PJ 




fc{z,r^) 


fR{z,r±) 









spin-isospin projection operator, and and are cou- 
pling strengths adjusted to the data. 

As seen in (1-84) and (1-89), for the commonly used 
pairing force Ijl33p . only two pairing densities come into 
play: the isovector density p and the isoscalar p-p spin 
density Sq. The corresponding isovector p-p potential 

U{r) is simply proportional to p while the isoscalar p-p 
field So(r) is proportional to the scalar product of the 
quasiparticle's spin & and Sq. Physically, p represents the 
density of S=0, neutron-neutron, proton-proton, and p-n 
pairs while the vector field Sq describes the spin distri- 
bution of 5=1 isoscalar p-n pairs. Indeed, when express- 
ing Eq. (|12bp directly in terms of the p-n pairs, one can 
see that Sqz contains the S—1, Ms—0 component while 
Sqx and SQy contain combinations of Ms=l- and Ms—~1 
pairs. The physical interpretation of the isoscalar p-p 
mean-field Hamiltonian 

^o(^) = So • o- oc Sq • o- (136) 

is the projection of the quasiparticle's spin on the spin of 
the p-n pairing field, its local magnitude determined by 
the HFB equations; hence, the SCSs are present in the 
problem. 

It is important to emphasize that the isoscalar density 
So contains all magnetic components of the S=l p-n pair- 
ing field. When studying an individual component sepa- 
rately, e .g., in the context of the so-called a — a ot a — a 
pairing |39j , one may arrive at erroneous conclusions that 
the presence of isoscalar pairing must be associated with 



breaking certain SCSs, such as axial symmetry or signa- 
ture. The usual argument, made originally in Ref. [39| 
and then repeated in the literature [40|, |41|, |42| , is that 
the individual components of the S=l pair field are not 
invariant under rotations. For instance, the pairing ten- 
sor KiM does not commute with signature TZa {a=x, y, z) 
0: 

7e-ViM„7ea = (-l)*^"«:iA./„, (137) 

and this has led to a conclusion that the isoscalar pairing 
must break signature. 

Let us consider axial and mirror symmetry as SCS. As 
seen in Table ITVl the isoscalar pairing density Sq vanishes 
only if the p-n symmetry is conserved. In the generalized 
pairing theory that mixes proton and neutron orbits, the 
solenoidal field Sq is nonzero. The lines of field So are 
schematically shown in the left panel of Fig.[TJ The right 
panel shows that while an individual vector Sq at a given 
point is not invariant with respect to symmetries S such 
as rotations around the third axis or signature TZ^-, the 
field So is perfectly covariant 

s^{r) = So{S+rS). (138) 

The scalar product (|136p is actually invariant in both 
cases shown in Fig. [TJj. It is interesting to note that 
for the geometry of Fig. [Tl the third component Soz as- 
sociated with the M=0 isoscalar pairing field vanishes. 
That is, the solenoidal pairing field is created by the two 
components with M—±l. Therefore, we conclude that 
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TABLE VII: Properties of local particle-hole rotationally symmetric (S0(3)-invariant) transition densities, depending on the 
conserved (C) or broken (B) space- inversion (P), proton- neutron (p-n), or time-reversal (T) symmetries. Generic real, imagi- 
nary, or complex functions of the radial variable r are denoted by fnir), fi{r), or fc{r), respectively. 



conserved (C) or broken (B) 



B 


B 


B 


B 


u 


(J 




(J 


g 










g 


n 


n 


B 


C 


B 


c 


B 


c 


B 


C 


fc (r) 


fair) 


fc{r) 




fc{r) 


fR[r) 


fcir) 


fR{r) 


fc (r) 










fc{r) 


fR{r) 








fc{r) 


fiir) 








fcir) 


fiir) 








fc(r) 


fR{r) 


fc{r) 




fc{r) 


fRir) 


fc(r) 


fR{r) 


fc{r) 


fR{r) 








fc{r) 


fR{r) 








fc{r) 


fi{r) 








fc{r) 


fiir) 








fc (r) 


fB{r) 


fcir) 


fR{r) 














fc (r) 


fR{r) 




















fc (r) 


fi{r) 




















fc{r)er 


fi{r)er 


fc{r)er 


fi{r)er 














fc(r)er 


fi{r)er 




















fc(r)er 


fR{r)er 




















fc(r)er 


fi{r)er 


fc{r)er 


fi{r)er 














fc(r)er 


fi{r)er 




















fc{r)er 


fR{r)er 




















fc{r)er 


fi{r)er 


fc{r)er 


fiir)er 














fc(r)er 


fi{r)er 




















fc{r)er 


.fi?(r)er 




















fc{r)er 


fiir)er 


fc{r)er 


fiir)er 


fcir)er 


fiir)er 


fcir)er 


fi{r)e 


fc{r)er 


fi{r)er 








fc{r)er 


fi{r)er 








fc{r)er 


fR{r)er 








fc{r)er 


fR{r)er 








fc{r)er 


fR{r)er 


fc{r)er 


/«(r)er 


fc{r)er 


fR{r)er 


fcir)er 


fR{r)e 


fc{r)er 










fc{r)er 


fR{r)er 








fc{r)er 


fiir)er 








fc{r)er 


fi{r)er 








fc{r)S 


/«WS 


fc{r)S 


fR{r)S 














fc(r)S 


fR{r)S 




















fc{r)S 


fiir)S 





















the assumption of axial symmetry, or signature, does not 
preclude the existence of isoscalar pairing. 

For the finite-range pairing forces and for the general 
Skyrme p-p functional, other pairing densities appear 
in addition to p and Sq. In the order of importance, 
the next crucial densities are the isoscalar vector cur- 
rent density jo^r) and the isovector spin-current tensor 

density J(r). Both can be generated by the momentum- 
dependent term in Eq. (|135[) and can be associated with 
the L=l pairing field. 

The isoscalar vector current density is solely responsi- 
ble for the isoscalar pairing field in the spherical, mirror 
symmetric case when time-reversal symmetry is broken. 
Indeed, according to Table |TT1 the L—0 density Sq van- 
ishes in this limit. It is only when the mirror symmetry 
is broken that sq becomes nonzero in spherical nuclei. 
Such a scenario could be an interesting possibility in very 
neutron-rich nuclei, in which, e.g., the l7rd5/2 and lvf5/2 
or l7rf7/2 and lvgj/2 orbitals could appear near the Fermi 

surface. The isovector tensor density J{r) is generally 



nonzero in spherical nuclei, and it has radial character. 

In the axial, parity-conserving case jo{f) is perpen- 
dicular to So- It is interesting to note that when the 
mirror symmetry is broken, the densities Sq and jo have 
all components (radial, azimuthal, and vertical) nonzero 
but their geometries will differ in general. 

In the absence of spin polarization, i.e., for time- 
reversal invariant systems, the isoscalar pairing field Sq is 
purely imaginary. In the presence of rotation, the time- 
reversal symmetry is internally broken. In this case, the 
pairing field is generally complex. Consequently, assum- 
ing the real Bogoliubov transformation and real pairing 
tensor 40] may limit the domain of self-consistent solu- 
tions. 

Another factor that may impact the generality of con- 
clusions of Ref. [i^l is the lack of the p-n symmetry- 
breaking on the HF level. Such an approximation does 
not seem to be justified as the self-consistent polarization 
between p-h and p-p channels is well known in the isovec- 
tor pairing case. Originally, the condition that the p-h 
density matrix preserves the p-n symmetry has been pro- 
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TABLE VIII: Properties of local axially and mirror symmetric (0^^(2)-invariant) particle- hole transition densities depending 
on the conserved (C) or broken (B) proton-neutron (p-n), or time-reversal (T) symmetries. The z-simplex (Sz) symmetry 
is conserved. The vector and pseudovector densities take the form of an expansion in a basis of some of the three vectors: 
r±, z and z x r±. The pseudotensor densities are linear combinations of the two symmetric pseudotensors: r± {z x r±) 

and z ^ {z X r±^). The expansion coefficients are real, imaginary, or complex functions of the two 0^^(2) scalars, 2^ and r']_. 
Generic real, imaginary, or complex coefficients are denoted by /_B(z^,r5^), fi{z^,r']_), or fc{z^,r'{^), respectively. 



G \ rm TTi o^" y\ t 
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> 
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fc{z\rl) 
fc{z\rl) 
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rr(t) 
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z X r± 
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posed in Rcf. [39] in the context of BCS calculations for 
N=Z nuclei where it was postulated that the expectation 
value of isospin in the quasiparticle vacuum is zero: 

(^f|f|*)==0. (139) 

We note that while the absence of the p-n mixing in the 
p-h sector automatically guarantees the condition (|139p 
for the N—Z nuclei in the absence of isospin-breaking 
interactions, the three constraints (|139p are in general 
not sufficient for all the t' — ~t matrix elements of p to 
vanish. 

In addition, the independent treatment of time- 
reversal and isospin symmetries as done in Ref. (39j is not 
justified. Indeed, as pointed out in Ref. [l^], the time- 
reversal and the isospin rotations do not commute. This 
implies that the relative phases between proton and neu- 



tron wave functions in a p-n broken quasiparticle state 
cannot be chosen arbitrarily. 

VII. CONCLUSIONS 

In this study, we investigated the symmetries of nucle- 
onic densities of the generalized nuclear DFT that allows 
for the arbitrary mixing of protons and neutrons. We 
considered the most important self-consistent symme- 
tries: spherical, axial, space-inversion, and mirror sym- 
metries. The main conclusions of our work can be sum- 
marized as follows: 

1. The local pairing densities po (isoscalar pairing den- 
sity) and s (isovector spin density) always vanish. 

2. One can always construct a phase convention for 
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TABLE IX: Properties of local axially symmetric (S0^(2)-invariant) particle-hole transition densities depending on the con- 
served (C) or broken (B) proton-neutron (p-n), or time-reversal (T) symmetries. The 2;-simplex (S^) symmetry is broken. 
The vector, pseudovector, or pseudotensor densities take the form of an expansion in the vector {e±,e^,ez) or tensor 
(P'^,^^^, S^, S^"*^) basis. The expansion coefficients are real, imaginary, or complex functions of cylindrical coordinates r± 
and z. Generic real, imaginary, or complex coefficients are denoted by fR{z,rj_), fi{z,r±), or fc{z,rj_), respectively. 
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which the local p-h densities are purely real [22| . In 
the presence of particular SCSs, some p-h densities 
vanish. 

3. In the absence of SCSs, the local p-p densities are 
complex. If time reversal is SCS, p-p densities be- 
come either purely real or purely imaginary. 

4. If p-n symmetry is SCS (no explicit p-n mixing), 
the k—1 and 2 isospin components of p-h densities 
and k—0 and 3 isospin components of p-p densities 
vanish. 

5. When 0(3) is SCS (spherical, mirror-symmetric 
case), the local pseudoscalar, pseudovector, and 
pseudotensor densities vanish. The only nonzero 
isoscalar-pairing density is the current density 
jo(r). All those densities can become nonzero in 
the SO (3) limit when the parity is broken. See Ta- 
bles H] and [TT] for a summary. 



6. When S0(2) and 5*2 are SCS (axial, mirror- 
symmetric case), pseudoscalar and pseudotensor 
densities vanish. If time reversal is SCS, p-p densi- 
ties become either purely real or purely imaginary. 
Properties of local axially symmetric (0^^(2)- 
invariant) densities are listed in Tables [1111 and [TVl 

7. When SO (2) is SCS (axial case), all densities 
are generally present. Properties of the S0(2)- 
invariant local densities are listed in Tables |V] and 

ED 

8. When space inversion, three signatures, and three 
simplexes are SCSs (D£^ group), all densities are 
generally present. Symmetry properties of p-h den- 
sities have been previously discussed in Refs. [13, 
[2^ . The analogous expressions for p-p densities are 
given in Sec. IIVCI 

9. The symmetry properties of the transition densi- 
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FIG. 1: (Color online) (a) Schematic illustration of the 
isoscalar vector field so in the case of conserved axial and 
mirror symmetries. The field is solenoidal, with vanishing 
third component, (b) Under rotation around the third (sym- 
metry) axis, the field at point rA is transformed to position 
rs. Likewise, under TZx, the field is transformed to rc- While 
neither of these operations leave the individual vector sqIta) 
invariant, the field as a whole does not change, i.e., it is co- 
variant. 



strong dynamical coupling between p-h and p-p chan- 
nels is expected. Consequently, in order to fully benefit 
from the p-n symmetry-breaking mechanism, p-n sym- 
metry should be broken already on the level of p-h mean 
field. This is not what has usually been done in existing 
calculations. 

The main results of this paper, summarized in Tables I- 
fX and in relations in Sec. lIVCl are symmetry properties 
of the local p-h and p-p densities that are building blocks 
of the generalized nuclear DFT formalism. These results 
can be useful when building a microscopic framework, 
rooted in the LDA, to describe various phenomena oc- 
curring in N'^Z nuclei. 
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(UNEDF SciDAC Collaboration), DE-FG02-96ER40963 
(University of Tennessee), DE-AC05-00OR22725 with 
UT-Battelle, LLC (Oak Ridge National Laboratory), and 
DE-FG05-87ER40361 (Joint Institute for Heavy Ion Re- 
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APPENDIX A: GENERALIZATION OF THE 
CAYLEY-HAMILTON THEOREM 



ties are the same as those of the local densities. 
The only notable difference is that the p-h tran- 
sition densities can be, in general, complex unless 
the time reversal is SCS and some additional re- 
strictions are present. 

10. The isoscalar pairing density Sq is the main build- 
ing block of the T=0 pairing field. In axial, 
reflection-symmetric nuclei, this field is solenoidal 
and gives rise to a paring potential that preserves 
signature. 

11. The second most important building block of the 
isoscalar pairing field is the current density Jq. It 
is solely responsible for the isoscalar pairing field 
in the spherical, mirror symmetric case when time- 
reversal symmetry is broken. For axial systems, 
this field is solenoidal. 

12. Breaking space inversion or mirror reflection may 
have profound consequences for the existence of 
isoscalar pairing as many isoscalar densities van- 
ishing in a parity-conserving limit can become 
nonzero. 

The symmetry properties discussed in this work pro- 
vide the necessary, but not sufficient, conditions for the 
presence of isoscalar pairing in nuclei. Whether such 
fields will appear or not depends, of course, on the actual 
form of the EDF and the values of coupling constants. 
In general, similar to the isovector pairing channel, a 



When analyzing symmetries, we often meet the prob- 
lem of how to construct a tensor quantity in terms of 
another tensor. Such a situation has been first encoun- 
tered in investigations of phenomenological constituent 
equations for macroscopic systems (see, e.g., [43|). The 
startingpoint in the analysis is the Cayley-Hamilton the- 
orem j44 . [45| by which an arbitrary function f{A) of a 
3x3 matrix A can be expressed as 

f{A)ab = Co{A)5ab + CM)Aab + ^2(^4)^2,, (Al) 

where the scalar functions Ci{A) = ^(01,02,03) {i—Q, 
1, or 2) depend only on three independent invariants of 
A: fli (trace of A), 02 (trace of A^), and 03 (determinant 
of A). In this case, matrix A is a rank-2 (reducible) 
tensor (A = 2) and the tensor field /(A) is also a rank-2 
(reducible) tensor [L = 2). Coupling of rank-2 tensors to 
L — 2 simply corresponds to multiplying matrices, and 
then Eq. (gl]) can be derived from the Taylor expansion 
of /(A) combined with the observation that every matrix 
obeys its own characteristic equation (any power A" (n > 
2) can be written as a linear combination of A°, A^, and 
A2). 

In applications presented in this study, we are inter- 
ested in spherical (irreducible) tensors and tensor fields 
corresponding to the rotation group 0(3), 

d(a/37) = e''-'^' e'^^^ e'"^- , ( A2) 

cf . Eq. (|43ap . In this case, xa/^ is called irreducible spher- 
ical tensor of rank A fiH] if in the rotated reference frame 
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it can be expressed as: 



where D^, ^{aP^) are the Wigner functions |46 
Now let $LM be a tensor field of xa^: 



*LM'= E ^M'Af(«/37)<I' 



LA/, 



(A3) 



(A4) 



A[=-L 

which is a function of 2A + 1 components of Xa/^j 

*LA/ = $LM (xam) , (A5) 

such that 



(A6) 



In other words, we are interested only in tensor fields 
being the isotropic functions of xa/j 143] ■ 

Condition (IA6|) is essential: it states that the rotated 
tensor field can be obtained by calculating the orig- 

inal tensor field ^lu at arguments which are the rotated 
tensor components x'a^' of X\^^■ It means that functions 
(jASp . apart from depending on XA/i, do not depend on 
any other tensor object (fixed material tensor). If they 
did, values of rotated tensor field could have been 

obtained by rotating the arguments xx^l ^'^'^ all mate- 
rial tensors simultaneously, but otherwise the rotation of 
arguments suffices. 

From Eq. (jAl[) we derive the general form of the 
quadrupole (L = 2) field $2Af b eing an isotropic func- 
tion of the quadrupole tensor X2 [43| : 

^'2Af(X2) = C'l([X2 X X2]o, [X2 X X2 X X2]o)X2Af 

+ C'2([X2 X X2]o, [X2 X X2 X X2]o)[X2 X X2]2Af- 

(A7) 

The notation [ x means the vector coupling to multi- 
polarity L. Symbols x and • for the vector and scalar 
products of vectors used previously in the paper are, up 
to coefficients, equivalent to [ x ]i and [ x ]o, respectively. 

To generalize Eq. (|A7p to arbitrary values of L and A, 
one should first establish a (finite) complete system of 
i\ irreducible elementary tensors eim characteristic for a 
given A. (Irreducible means that none of them can be 
expressed rationally and integrally in terms of the oth- 
ers.) The highest components (m = I) of the elementary 
tensors are called elementary factors. The elementary 
tensors are constructed by successive vector couplings of 
n xa's to different intermediate Vs (meaning that I < nX): 



- An) 



(H) = [XA x • • • X xa; 



[c] 



(A8) 



e^^^ is itself an elementary tensor. In the particular case 
of Z = 0, all independent scalars are elementary factors. 
Completeness of the system of elementary factors does 
not exclude the relations (syzygies) between them. The 
syzygies can be written in the form: 



^3 i^li ' ■ 



(A9) 



for j ~ I, . . . , jx, where Sj are rational integral functions. 
Again, the number of independent syzygies is finite and 
depends on A. 

Having determined the elementary tensors, we align 
them (i.e., couple to the maximal multipolarity) to get a 
tensor of rank L: 



e, X 



\L=l+V 



(AlO) 



Alignment of the elementary tensors means multiplica- 
tion for the elementary factors. Because of the exis- 
tence of syzygies ()A9|) . some aligned tensors can be ex- 
pressed in terms of others. Using syzygies, we find a fi- 
nite number k\L of independent aligned tensors (y9|j^^(e), 
(k — l,...,k\L) called fundamental tensors. It turns 
out that an arbitrary tensor field of rank L being the 
isotropic function of tensor xa can always be presented 
in the form: 



LM 



(xa) = ^Cfe(eo)v'iA/(£), (All) 



fe=i 



where argument Eq of Ck stands for all the independent 
scalars. Scalar functions can be, in general, arbi- 
trary. However, the form of some Ck can be restricted. 
For instance, high powers of some scalars do not ap- 
pear. This is because that a s^zygy can make an expres- 
sion (£oo)"(4';')" • • • for some i and i' , . . . ,i" , 
and n and n', . . . ,n" dependent on other elementary fac- 
tors. We refer to Eq. (jAll|) as the Generalized Cay ley- 
Hamilton (GCH) theorem. 

An explanation of the procedure presented above can 
be traced back to the theory of covariants of algebraic 
forms given by Dickson [48], which, however, uses quite 
different notations than those used here. Let us sketch 
the main points of the theory 44] . Let the 2A-ic (of order 
2A) algebraic form, binary in variables a;i, X2^ be given 
by: 



^2A(a;i,a;2;xA-A, ■ • ■ ,xaa) 

2A 

A — /.J 



E f xl^,V^^A.x^%^-^ (A12) 



/i— — a 

with xa^ (m = ^'^1 ■ ■ ■ being the set of coefficients. 
Replacing the variables by a linear non-singular transfor- 
mation 



for i = 1,...,^^, where symbol [c], redundant in most 
cases, stands for a specific coupling scheme. Tensor xa = 



E ^ikVk 



(A13) 



k=l 
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for i = 1,2, in the form of Eq. (IA12p . one obtains: 

F2x{xi,x2;x\) = G2\{yi,y2;x\) 

(A14) 

where -0^/^ is a new set of coefficients. The 2Z-ic form, 

77l = -l ^ ^ 

(A15) 

such that 

H^;\yi,y2;M = (det(A,fc))'"if^r^(2:i, 2:2; xa) (A16) 

we call a (homogeneous) covariant of weight w of F. Co- 
efficients are homogeneous polynomials of order n 
(called the degree of the covariant) such that Xn = l + w. 

When / = 0, i?^"^ = /i^q^ is an invariant of F. The poly- 
nomial in front of the highest power of xi in Eq. (jA15[) . 
/ij"^ is called a semi-invariant of F. The theory of covari- 
ants of algebraic forms shows that an arbitrary covariant 
can be expressed rationally and integrally in terms of a fi- 
nite, irreducibly complete set of covariants, which can be 
related rationally and integrally to each other by a finite 
system of independent syzygies (Gordan-Hilbert Finite- 
ness Theorem, see (43|). The number of basis covariants 
for a few lowest A's is listed by Olver [ioj . 

What do the covariants of algebraic forms have in com- 
mon with the tensor fields as functions of tensors? ft 



turns out that semi- invariant /i;" (xa) forms the highest 
projection of a tensor of rank I dependent on the tensor 
of rank A (the heaviest state of an irreducible representa- 
tion of SO (3) embedded in an irreducible representation 
of SU(2A + 1)). The remaining polynomials in Eq. (|A15[) 

^Im (xa) are other components of the same tensor. Proof 
of this statement is outlined in Ref . [sO] . 

Constructive proof of Eq. (|A11[) can be performed by 
an explicit construction of a basis in the space of func- 
tions of x\ and demonstration that it has the structure 
of Eq. (jAlip . In case of A = 1, we know that an arbitrary 
tensor field of rank L as a function of the position vector 
r, takes the form: 

*LAf(r) = CL{r)YLM (^) , (A17) 

where = r r and Y^m is the spherical harmonic. From 
Eq. (|A17p we see immediately that in case of A = 1 , there 
are two elementary tensors {ii = 2), namely vector r and 
scalar r^, and no syzygy (ji — 0). For every given L there 
is only one fundamental tensor (fci^ = 1) of the form [4^ : 

V'lmI'^) = [r X r X ■ ■ ■ X vJ lm 

L 

(X r^^LM f-) . (A18) 
In particular, Eq. (|A17P tells that an arbitrary isotropic 
vector field takes the form: 

*(r) = C(r)r. (A19) 

Other examples are constructions of the oscillator bases 
in cases of A = 2 (see Ref. [131 and references quoted 
therein) and A = 3 [50]. 
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